Astronomy & Astrophysics manuscript no. sandin 


© ESO 2008 


May 23, 2008 





oo 

o 

O 

(N 

>> 



m 
<N 

6 



(N 
> 

00 

m 

(N 
O 
OO 
O 



- 

X 



Spatially resolved spectroscopy of planetary nebulae 

and their halos * 

I. Five galactic disk objects 

Christer Sandin, Detlef Schonberner, Martin M. Roth, Mattiiias Steffen, Petra Bohm, and Ana Monreal-Ibero 

Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-144 82 Potsdam, Germany 
Received 25 February 2008 / Accepted 20 May 2008 

ABSTRACT 

Strong mass loss off stars at the tip of the asymptotic giant branch (AGB) profoundly alfects properties of these stars and their 
surroundings, including the subsequent planetary nebula (PN) stage. With this study we wanted to determine physical properties of 
mass loss by studying weakly emitting halos, focusing on objects in the galactic disk. Halos surround the, up to several thousand times, 
brighter central regions of PNe. Young halos, specifically, still contain information of the preceeding final mass loss stage on the AGB. 
In the observations we used the method of integral field spectroscopy with the PMAS instrument. This is the first committed study of 
halos of PNe that uses this technique. We improved our data analysis by a number of steps. In a study of the influence of scattered 
light we found that a moderate fraction of intensities in the inner halo originate in adjacent regions. As we combine line intensities 
of distant wavelengths, and because radial intensity gradients are steep, we corrected for effects of differential atmospheric refraction. 
In order to increase the signal-to-noise of weak emission lines we introduced a dedicated method to bin spectra of individual spatial 
elements. We also developed a general technique to part the temperature-sensitive oxygen line [Oiii] /14363 from the adjacent telluric 
mercury line Hg/i4358 - without using separate sky exposures. By these steps we avoided introducing errors of several thousand 
Kelvin to our temperature measurements in the halo. For IC 3568 we detected a halo. For M 2-2 we found a halo radius that is 2.5 
times larger than reported earlier. We derived radially densely sampled temperature and density structures for four nebulae, which all 
extend from the central regions and out into the halo. NGC 7662, IC 3568, and NGC 6826 show steep radially increasing temperatures 
and a hot halo, indicating that the gas in the halo is not in thermal equilibrium. M 2-2 shows a larger temperature in the central region 
and an otherwise constant value. From the density structures we made estimates of core and halo masses and - for the first time 
reliable - mass loss rates at the tip of the AGB. All four objects show inwards radially increasing mass loss rate structures, which 
represent a rise by a factor of about 4-7, during the final mass loss phase, that covers a time period of approximately 10"* years. Within 
a factor of two, the average of the maximum mass loss rates, which are distance dependent, is Mniax - lO""* MQyr"' . 

Key words. Methods: data analysis. Techniques: spectroscopic, stars: mass-loss, planetary nebulae: individual (IC3568, M2-2, 
NGC 3587, NGC 6826, NGC 7662) 



^ 1. Introduction mass loss stage is short, and selection therefore works against its 

detection, making detailed observational studies difficult. 

Stars on the thermally-pulsing asymptotic giant branch (AGB) 
H are in the final stages of evolution before they turn into cooling Planetary nebulae (PNe) are in transit between the tip of the 
■ - white dwarfs. During this phase these stars lose nearly all of the ^GB and the white dwarf domam of the Hertzsprung-Russell 
mass outside their core in an increasingly strong stellar wind, diagram. These extended objects occur in rather large numbers 
The matter of such winds is enriched with nuclearly synthesized ^an readily be observed in the visible wavelength range, 
elements, which have been dredged up from the interior during where a significant fraction of the luminosity is emitted in a 
the evolution; thereby contributing to the galactic chemical evo- f^w emission hnes. The morphological variation of PNe is large 
lution. Models of strong (dust-driven) stellar winds are described (^ee, e.g., Balick & Frank 2002; Sahai et al. 2007, and refer- 
by, for example Bowen (1988), Fleischer et al. (1992), Winters ^nces therein). Of those nebulae with a spherical or elliptical 
et al. (2000), Hofner et al. (2003), and Sandin (2003, 2008). ^^ape the central parts consist of a relatively bright rim with a 
Two recent overviews of mass loss from cool stars are given by f^i"ter attached shell; both are built-up by an interplay between 
Willson (2000) and Hofner (2005). Although there is progress photo-ionization and wind interaction (see, for example, Frank 
in understanding both stellar evolution and mass loss theoreti- ^t al. 1990; Marten & Schonberner 1991; Frank 1994; Mellema 
cally, observational details of, in particulai-, the last phase at the 1994; Perinotto et al. 2004). Often a faint, and in many cases rel- 
tip of the AGB have remained obscure. The lifetime of the final ^tively large, halo surrounds the central PN region; see the PN 

halo atlas of Corradi et al. (2003, hereafter CSSP03) for exam- 



Send offprint requests to: Christer Sandin, e-mail: CSandin@aip . de Ples- Such halos are with the help of models identified as the 

* Based on observations collected at the Centro Astronomico lomzed AGB wind (Steffen & Schonberner 2003; Villaver et al. 

Hispano Aleman (CAHA), operated jointly by the Max-Planck Institut 2002a,b). Halos of younger nebulae are less affected by PN shap- 

fur Astronomie and the Instituto de Astrofi'sica de Andalucia (CSIC). ing processes than the central parts are, and therefore contain in- 

** Figures 24-29 are only available in electronic form via formation about the preceeding mass loss episode; this is not the 

http://www.edpsciences.org. case with so-called recombination halos of older nebulae (see. 
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e.g., Tylenda 1986; Corradi et al. 2000). Due to very low surface 
brightnesses little is so far known about physical properties of 
halos. The purpose of this study is to improve our knowledge by 
studying halos of a number of PNe in detail. 

Integral field spectroscopy (IFS) is a suitable tool for per- 
forming spectroscopy on regions spanning very weak to very 
strong surface brightnesses, such as PNe with halos. Through 
a generally large field of view IFS permits a co-adding (bin- 
ning) of any number of spectra of individual spatial elements 
on an integral field unit (IFU), which covers a defined area 
on the sky. In this study we used the Potsdam Multi-Aperture 
Spectrophotometer (PMAS, Roth et al. 2005) in order to per- 
form a plasma diagnostic study of selected PNe, including both 
central parts and halos. Our goal was to make measurements of 
electron temperatures and densities, in order to provide spatially 
resolved radial structures of these properties. While a density 
allows a determination of mass loss rates, a knowledge of the 
electron temperature is important when calculating abundances. 

We discern between objects in terms of the metalicity, and 
use four groups of objects: galactic disk, halo, Magellanic Cloud 
(MC), and other extragalactic PNe. The five objects we selected 
for this study are all galactic disk nebulae with circular outer re- 
gions. NGC 7662 and NGC 6826 are both well studied objects of 
considerable size with a well-confined halo. IC 3568 is a closely 
circular nebula at intermediate distance, that, so far, is reported 
not to have a halo. M 2-2 is a small and distant round nebula, that 
is not often studied. Although a present halo is reported in the Ut- 
erature, the center-to-halo flux ratios are too low (2-5, Guerrero 
et al. 1998, hereafter GVM98) to verify the existence of a weak 
halo (these ratios must be > 100, CSSP03). NGC 3587, finally, 
is a large evolved nebula with a potential recombination halo. 
The entire radial extent of the halo was, furthermore, not cov- 
ered for all objects, but the shell-halo transition region and inner 
halo were covered in every case. Priority was given to cover the 
blue wavelength range, including B aimer lines and the temper- 
ature sensitive oxygen lines. Some observations were addition- 
ally performed in the red wavelength range, in order to cover Ha; 
and the density sensitive sulfur lines, where present. It should be 
noted that only a few emission lines are visible in the halo, where 
densities normally are too low to be determined using Une ratios. 

The IFS method allows a determination of physical proper- 
ties in two dimensions, but studies presenting two-dimensional 
maps of PNe are not new. Temperature maps of various neb- 
ulae, using different methods of observations are presented by 
Reay & Worswick (1982), Jacoby et al. (1987, hereafter JQA87), 
Mavromatakis et al. (2001), Hyung et al. (2001), Rubin et al. 
(2002), and Wesson & Liu (2004). Liu et al. (2004) and Monteiro 
et al. (2004, 2005) derive densities and temperatures for two 
nebulae using a long slit scanning method. Tsamis et al. (2008, 
also see Tsamis et al. 2007) measure temperatures, densities and 
abundances for three nebulae using the FLAMES instrument and 
the Argus IFU. All these studies focus on the central, brighter, 
parts of nebulae. In the only study that is similar to ours, so far, 
Martm-Gorddn et al. (2007), using PMAS with the PPAK IFU, 
calculate a temperature map of a knot in the halo of NGC 6543. 

In the following we give the details of the observations in 
Sect. 2, followed by a description of the performed data reduc- 
tion in Sect. 3, where we also describe how we correct for dif- 
ferential atmospheric refraction (DAR). Before we present the 
results in Sect. 5, we treat important issues of the data analysis 
in Sect. 4 - including a discussion on scattered light and descrip- 
tions of new dedicated methods to deal with sky subtraction and 
binning of spatial elements. Our interpretation of the physical 



results are then found in Sect. 6. The paper is closed with con- 
clusions in Sect. 7. 

2. Observations 

On five difi^erent occasions the five objects studied here were ob- 
served with the PMAS instrument, that is attached to the 3.5m 
telescope at Calar Alto. We used the lens array (LARR) and 
PPAK IFUs (Roth et al. 2005; Kelz et al. 2006). Instrumental 
configurations, weather conditions, and details on the exposures 
taken are given in the observational journal in Table 1. In ad- 
dition to the science exposures continuum and arc lamp flat- 
fields were taken, as well as spectrophotometric standard-star 
exposures. In order to correct for varying fiber-to-fiber trans- 
mission sky flat-fields were taken for all targets, but NGC 6826, 
where a normaUzed continuum lamp flat-field was used instead. 
Note that in most cases no separate observations of the sky - 
completely void of halo emission lines - were made; instead 
an alternative method was adopted where a set of telluric lines 
are fitted together with the object lines, cf. Sect. 4.2. Data of 
NGC 3587 from 2004 were sampled using the beam switch tech- 
nique (BSW, cf. Sect. 4.2). This data, however, did not permit an 
accurate detection of the weak auroral line [O iii] /14363, and the 
object was therefore observed again in 2006, then without using 
the BSW technique. 

The LARR IFU holds 16x16 separate fibers, where each fiber 
represents a spatial element, or a so-called spaxel, on the sky. In 
the r.'O (0'.'5) sampling mode every LARR pointing (IFU) covers 
an area of 16"xl6" = 256"^ (8"x8" = 64"^) on the sky. In 
comparison the PPAK IFU holds 331 fibers, which are arranged 
in a hexagonal grid (see fig. 5 in Kelz et al. 2006), and cover 
an area on the sky of about l.'O^. Any number of spaxels can be 
co-added to create a final spectrum (cf. Sect. 4.3). 

We show images of the objects we studied in Fig. 1, where 
the locations of individual IFU pointings are indicated. These 
images can, with the exception of M 2-2 and IC 3568, be com- 
pared with the respective image in CSSP03 (the image of 
NGC 6826 is taken from this reference). 

Wavelength regions and gratings were chosen to cover emis- 
sion lines crucial when making a plasma diagnostic. All ob- 
jects were observed in the wavelength region covering the oxy- 
gen lines [Oiii] /14363, /14959 and (in most cases) A5007, and 
the Balmer lines Hy and HyS. Some of the other emission Unes 
in this blue part of the spectrum, observed simultaneously, 
were [NeiiiU3869, H<5, He 1/14471, Heiii4686, and the den- 
sity sensitive doublet [Ariv] AAAIW, 4740. Note, however, that 
[Ariv]/I4711 is blended by He 1/14713, which can make the 
use of this line ratio uncertain (see sect, n in StanghelUni & 
Kaler 1989, hereafter SK89, and sect. 2 in Wang et al. 2004). 
Moreover, [Ne iii] /13968 blended with He was detected in ev- 
ery case when [Nem]/i3869 was observed, but was not con- 
sidered further. A resolved density sensitive oxygen doublet 
[O ii] /1/13726, 3729, requiring a separate setting with a high res- 
olution grating, was only observed for two of the five objects. In 
the red part of the spectrum the goal was to cover Ho-, and the ni- 
trogen fines [N n] /1/16548, 6583, along with the density sensitive 
sulfur doublet [S n] XX61\1, 6731. 

3. Data reduction 

The tools we used in the data reduction, P3d_online and 
Ppak_online, are parts of the PMAS P3d pipeline (Becker 2002, 
Roth et al. 2005). We first describe the data reduction method for 
observations carried out with the LARR IFU. 
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Table 1. Journal of observations 



Object 


Date 


Grating/ 


Wavelength 


IFU^ 


Weather conditions 


Exposures'^ 




KJUai^L VCU 


- [A/pixel] - 


[A] 




Sky Quality 


o ■ h 

Seemg" 


[seconds (s) or minutes (m)] 


NCjC 3587 


20U4.02.iz 


^71 r\n /I an 1 or\(\ 

V300/ 1.67/ 800 


icon nr\r\c\ 

3590-7000 


yir. BSW 


photometric 


1'.'1-1'.'4 


9Am:1.20 9r)m:l.05 ^j^m : 1.06,1.10 
^S4W::a' ^^'llOW::b' ^"l26W::c 




ZUU4.UZ.lJ 


V IZUU/U. Jo/ZoVU 


jo4U-^4UU 


i//riBSW 
"LARR 


photometric 


1''5-Z'0 


'^rjmil.lS •ir\m:1.10 '^rim:1.05 
'^"94W::a' ^"llOW::b' ^"l26W::c 




ZUUo.U4.Zo 


V 1ZUU/U.jD/jU8U 


OZIU-OO/U 


1"A 

i- Ularr 


not optimal 


r.'i-r/2 

rv'Q v'\ 
U. y— 1 . 1 


9f^m : 1 .05, 1 .05 orim : 1 .06, 1 .08 
^"74W::1 '^"90W::2 
^rim: 1.10,1.12,1.17 ^rim : 1.21,1.29 
"l06W::3 ' ^"122W::4 




2006.04.28 


V1200/0.35/3230 


4300-5040 


1^ • U LARR 


very good 


l'.'5-2'.'0 


^rjiTi: 1.07,1.09,1.13 -jnm: 1.16,1.23 oAm:1.37 
^"90W::2 ' "l06W::3 ' -^"l22W::4 


IC3568 


2004.02.16 


V300/ 1.67/ 800 


3590-7000 


y/ABSW 

^•"larr 


photometric 


FT 


cm: 1.42 Trim: 1.42 
-'c::0 '^"c::l 


NGC7662 


2005.09.09 


V600/ 0.81/1340 


3490-5150 


I'-'Olarr 


not optimal 


1'.'0-1'.'7 


ICS: 1.32 'iAAs:1.28 orym: 1.11,1.17 
'^-'c::l ' ^""l6W::2' -'"32W::3 


NGC6826 


2006.10.28 


V1200/0.35/3210 


4270-5010 


2'.' 68 ppAK 


very good 


1'.'6 


C!ris:1.06 cm: 1.08 QAm:1.09 
l :: I ' -^c:: 1 ' ■'"30SH::2 


M2-2 


2007.10.04 


V600/ 0.81/1360 


3550-5200 


0'.'5 larr 


photometric 


1'.'3-1'.'5 


,f.m: 1.16,1.16.1.16 . j-m : 1 . 1 1,1.08 4j-m:l.07 
-'"c::l '^-^8W::2 '^-^16W::3 




2007.10.05 


R600/ 0.78/1980 


5370-6960 


^.'5 LARR 


photometric 


1'.'0-1'.'3 


iQm: 1.07,1.06,1.06 ^cm:1.06 4cm: 1.08 
^'^c::l ' ^-'8W::2 ' ^-'l6W::3 



' The number indicates the size of each spatial element. The subscript indicates what IFU was used, and a superscript indicates whether beam switching (BSW) was used; cf. Sect. 4.2. 
^ The seeing is estimated from guide star measurements in the R-band. 

The superscript index (j: : x,x, ...) indicates minutes (m) or seconds (s), followed by the airmass of the respective exposure of a particular pointing, the subscript index (x :: x, ...) indicates 
the location (specified in arcseconds) of the center of the IFU relative to the position of the central star (c), followed by the tile number/letter used in Fig.l, 




Fig. 1. Object images illustrating where observations were done. The five panels show: a) NGC 7662, b) IC 3568, c) M 2-2, d) NGC 3587, and e) 
NGC 6826. The images in panels a-d were obtained with the PMAS Acquisition and Guiding (A&G) camera, and are reproduced in the same 
scale. The image of NGC 6826 in panel e was observed with JKT and is talcen from the image atlas of CSSP03. The bright disk visible in panels 
a and d corresponds to the diaphragm of the A&G pick-off mirror. The numbered squares/hexagons outline different positions of the PMAS 
IFUs across the respective FN shell and the halo. In panel c the inset shows a corresponding A&G R-filter image emphasizing the ansae in the 
southernmost part of the shell; cf. Sect. 5.3. In panel d the tiles drawn with dashed lines were observed on 12.02.2004 and 15.02.2004, and the 
tiles drawn with soUd lines on 26.04.2006 and 28.04.2006; cf. Table 1. 
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At first the bias level was subtracted and cosmic-ray hits 
removed. Second, a trace mask was generated from an inter- 
nal continuum calibration lamp exposure, identifying the loca- 
tion of each spectrum on the CCD along the direction of cross- 
dispersion. Third, a dispersion mask for wavelength calibration 
was created using an arc-exposure. In order to minimize effects 
due to a significant flexure in the instrument all continuum and 
arc lamp exposures were taken within two hours of the respec- 
tive science exposure. Fourth, a correction to the fiber-to-fiber 
sensitivity variations was applied by dividing with an extracted 
and normalized sky or continuum lamp flat-field exposure. In 
this process the data was changed from a CCD-based format 
to a row-stacked-spectra (RSS) format. Fifth, the data was cor- 
rected for differential atmospheric refraction (DAR), using the 
procedure outhned below. As a final step flux caUbration was 
performed in iraf using standard- star exposures. 

The reduction procedure is slightly different with PPAK data. 
Apart from a different number of spaxels and a different spaxel 
geometry, special care was taken to compensate effects of instru- 
ment flexure. Additional calibration fibers included in the PPAK 
fiber bundle were illuminated with a Thorium/ Argon lamp, pro- 
ducing single spots on the CCD. The locations of these spots 
shift with flexure in the same way as the measured spectra and 
can therefore be used to calculate shifts between calibration and 
science exposures. After fixing the so-called shift boxes the steps 
of the data reduction are the same as for the LARR IFU, i.e., 
make a bias, generate a trace mask from a continuum (halogen) 
lamp exposure, generate a dispersion mask from an arc-lamp ex- 
posure (mainly He-nHg), and generate a flat-field either from 
dome flat or continuum lamp exposures; we used a continuum 
lamp exposure. Afterwards the spectra were extracted from the 
science exposures by applying the prepared masks/frames with 
respect to the calculated shift boxes. The result is stored in the 
RSS format. 

In several cases, and in particular where there was only one 
science frame, cosmic -ray hits were removed using the L.A. 
Cosmic routine (van Dokkum 2001). Since this routine is adapted 
for images and not spectra, the input parameters had to be chosen 
carefully to ensure that ideally no, or as few spectral lines as pos- 
sible, were removed together with the cosmic -ray hits. Typical 
values used with the routine were: fn^ - 30-50 and (Tdip =4-12; 
the most careful (large) numbers were required with the standard 
star raw data. Any remaining hits at wavelengths close to emis- 
sion lines were then manually removed using splot in iraf on the 
final RSS-data - before correcting for DAR or flux calibrating. 

In order to capture all of the standard star flux, we summed as 
many spaxels as possible. These suimned spectra were then sky 
subtracted using an averaged sky spectrum, before creating the 
respective sensitivity function. This is not a problem since these 
spectra typically show few emission lines. Unlike the square 
lenslets, which are not separated on the LARR IFU, the fibers 
on the PPAK IFU are circular, of 2'.'68 diameter, and are sepa- 
rated by 3'.'6. Making an absolute flux calibration in this case, 
the total flux of the standard star has to be corrected for light 
lost between the fibers. We made such a correction by first fit- 
ting the intensity and position of a two-dimensional Gaussian to 
the standard star exposure data. Thereafter the fraction of light 
illuminating the area between fibers was estimated, taking the 
shape of the Gaussian curve into account. 

Several of the objects were observed at such a high air- 
mass that effects of DAR become visible - i.e., a wavelength- 
dependent spatial shifting of spectral images on the IFU. Other 
factors, which also determine the extent of effects of DAR are 
the geometrical size of the spaxels, and the presence of strong 



intensity gradients. In effect, the smaller the spaxels the more of 
the flux of different wavelengths end up on separate spaxels. It 
is important to correct for DAR when comparing lines of differ- 
ent wavelength. In order to do the correction, and (re)place flux 
from the same spatial region in the same spaxel, the approach of 
Filippenko (1982) is used to calculate an image offset vector for 
each wavelength bin. Each image, at every wavelength, is there- 
after shifted with respect to a pre-defined reference wavelength 
(that mostly is set to A^ef = 5050 A), using a (fractional) bilinear 
interpolation procedure. Doing this the intensity in each spaxel 
becomes a function of four spaxels - introducing an unavoidable 
smoothing, which additionally makes the intensity error estima- 
tion more challenging. In contrast to our approach Arribas et al. 
(1999) use shifting positions of an intensity maximum at dis- 
crete wavelengths in the data to calculate the offset vector. Their 
approach is difficult to apply when a maximum is located to the 
edge of the IFU, which for example is the case when the IFU is 
placed in the halo region where intensities decrease monotoni- 
cally. 

An additional consequence of the DAR correction procedure 
is that those strips of spaxels on the two edges of a (rectangu- 
lar) IFU have to be masked, which for some wavelength bins 
depend on flux from regions outside the IFU. The width of the 
strips on each side, and on which side they have to be masked, 
depend on the airmass during the observations, the wavelength, 
and the location on the sky; see Figs. 4 and 26 for examples. 
Note that no correction for DAR was made for NGC 6826, that 
was observed with PPAK at a rather low airmass. The large and 
separated spaxels of this IFU make the usefulness of the interpo- 
lation, as described above, doubtful. 



4. Data analysis - issues critical to results 

Before presenting final results in Sect. 5 we next address four 
topics that all have an important influence on properties of the 
outcome. We first discuss the influence of scattered light, esti- 
mated with a measured point spread function, in Sect. 4. 1 . Then 
we present a method for dealing with telluric (i.e., sky emission) 
Unes in Sect. 4.2. This is foUowed by the presentation of a dedi- 
cated method for binning spectra to increase the signal-to-noise 
in Sect. 4.3, and a discussion of dereddening in Sect. 4.4. 

4.1. Estimating tlie influence of scattered liglit 

In the halo of a planetary nebula the surface brightness in any 
line can be more than a thousand times weaker than in the rim, 
where most flux is emitted. An instrumental point spread func- 
tion (PSF) of a single source can be measured out to very large 
radial distances, where the relative intensity is 10^ times as weak 
as in the center (King 1971). In their work on halos of PNe, 
Middlemass et al. (1989a) and CSSP03 (sect. 4.1) study effects 
of scattered light. A qualitative study is motivated also for us 
to see where emission in the halo region originates, because all 
instruments scatter light to some degree. 

In order to determine the instrumental PSF we first observed 
the bright star Vega (orLyrae) in the red part of the spectrum, 
across three spatially offset tiles (with consecutively increasing 
exposure times), using the LARR IFU in the 0'.'5 setting. The 
maximum distance from the central star, where we measured an 
intensity, was thereby r - 20". We extracted a radial slice, one 
spaxel wide (O'.'S), covering all three tiles at the wavelength of 
Ha, and then used it to create an axi-symmetric two-dimensional 
PSF kernel of 40" diameter. The one-dimensional PSF is illus- 
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>, Fig. 2. Estimating effects of scattered light in 
M the observed data. Panel b shows the instru- 
B mental point spread function (PSF) measured 
-2^ for the bright star Vega, out to r = 20" away 
J from the PSF center. The ordinate is logarith- 
-3 £ mic. The three profiles that are vertically off- 
o set, by -1 dex incrementally, show the PSF trun- 
^ o) cated at: rpsF = 8" (dashed line), rpsF = 4" (gray 
— line), and rpsp = 2" (light gray line). The two 
tile boundaries are indicated with vertical dot- 
ted lines. The observed intensity structure of 
H/3, that is convolved with the PSF, is illustrated 
in panel a (solid black line), again with a log- 
arithmic ordinate. The observed normalized in- 
tensity structure is indicated with crosses and 
the light gray line. For a study of the appear- 
ance of a non-existent halo the structure was 
truncated at r= 15" (shown by the gray line that 
is vertically offset by -0.5 dex), and then con- 
volved with the PSF. The resulting structure is 
shown with a dash-dotted line. Convolved-to- 
observed ratios of the same two structures are 
shown in panel c. Panel d, finally, shows how 
the convolved structure plot depends on the ra- 
dial extent of the PSF. The thin black line uses 
the full PSF (as in panel c). For the remaining 
three lines we used a truncated PSF - the col- 
ors of these lines are the same as in panel b. 
The figures show that in the inner halo the con- 
_ 90 tribution of scattered light mostly comes from 

adjacent regions (gray lines; d). For further de- 
tails see Sect. 4.1. 
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trated in Fig. 2b; the measured full width at half maximum of the 
PSF is about l'.'4-r.'8. Note the jump in intensity at the transition 
between the two inner tiles - the amount of scattered light de- 
creases by about 0.6 dex when Vega is outside the field of view, 
i.e. outside the lens array and the diaphragm of the A&G camera 
pick-off mirror, cf. Fig. 1. In this context it should be mentioned 
that the most accurate PSF - which would extend across sev- 
eral tiles - ought to be determined for each spaxel individually, 
otherwise tile-to-tile jumps are placed incorrectly. Our aim was, 
however, to make a qualitative estimate of the influence of the 
PSF and we did not consider such a detailed treatment. 

Next we created a two-dimensional intensity structure, of 
80" diameter, using the radial surface brightness of NGC 7662 
in Hy6 (cf. Sect. 5.1.1; the seeing was about the same as dur- 
ing the night the PSF was observed). This structure was then 
convolved with the PSF kernel in order to achieve a first order 
estimate of effects of scattered light. Note that in reality the ob- 
served intensity structure must be understood as the convolution 
of the true surface brightness distribution with the same PSF. We 
show a radial slice of the intensity before and after convolution 
in Fig. 2a. The convolved-to-observed intensity ratio is shown in 
Fig. 2c, indicating the difference more clearly. While differences 
are small inwards of the shell-halo transition (i.e. for r $ 15"), 
there is a contribution of up to 30% to the intensity in the outer 
region from the central regions of the nebula (also see below). 

As a test, to see what the intensity structure of a non-existent 
halo would appear like we then truncated the two-dimensional 
observed intensity structure at r = 15", before convolving with 
the PSF kernel. The result, that is equivalent to a PN without 
a halo, is illustrated with the dash-dotted lines in Figs. 2a,c. 
Compared to the solid line the intensity decreases sharply for 



15"<r< 18", is ^ 15% of its strength, and decreases for r> 18". 
This demonstrates that the observed halo is real, because the ob- 
served data does not show such a steep decline. Likewise this 
test also shows that a fraction of the emitted intensity in the halo 
originates in the central regions; this fraction is the largest in the 
inner halo. 

In order to study to what extent different parts of the PSF 
contribute to the emitted intensity in the halo we made three tests 
where the PSF is truncated at increasing radii. Figures 2b and d 
show the result when such a truncation was made at rpsp = 2"0, 
'"PSF - 4'.'0, and rpsF - 8'.'0. Evidently local regions contribute 
more to the emitted intensity in the halo for 15"<r<18", where 
about two thirds of the scattered light comes from the inner 4" 
of the PSF; and the rest from the remaining PSF. In the outer 
region, for 18" < r < 30", it is in contrast the outer part of the 
PSF, beyond rpsp — 8", that contributes the most; the intensity 
is there ^ 10% larger because of scattered light from the outer 
wing of the PSF. 

Radial density structures, which are calculated using the 
Abel transform on intensities (cf. Sect. 5), are affected by scat- 
tered light. Estimating the influence we first calculated an ap- 
proximative deconvolved two-dimensional intensity structure of 
H/3 - instead of using a forwards convolution; by matching the 
convolution of the deconvolved structure and the PSF with the 
observed intensity structure we could determine it by iteration. 
The resulting density structure, that is based on this deconvolved 
structure, was found to be up to about 30 per cent less dense than 
the density structure, that we derived from the observed intensity 
structure (in the inner halo). This error estimate for the density 
is directly comparable to the amount of scattered light that we 
found in the intensity seen in Fig. 2c. 
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In comparison to these tests it is more difficult to see what 
the influence is on properties which depend on line ratios. We 
therefore also convolved the radial surface brightness structures 
of oxygen [O iii] /14363 and /15007, and found that the resulting 
temperature structure is lower by ^ 6%, compared to the mea- 
sured temperatures. The largest difference is seen in the halo. 
This contribution to the error of the temperature is of a similar 
magnitude as that due to a badly known extinction coefficient in 
the halo (cf. Sect. 5.1.2). 

Consequently, because of scattered light the observed sur- 
face brightness of any line is about 0-30% brighter than the 
true brightness, and the deviation is the largest in the inner halo, 
where gradients are the steepest. For line ratios of a particular 
region - and quantities which depend on these ratios, such as 
the electron temperature - this implies that actual emission is 
mixed with emission from regions of different physical prop- 
erties. True extrema in these regions are thus even more pro- 
nounced than what is observed. Since the differences are rea- 
sonably small, however, we did not correct our measurements 
for scattered light, but note that such a component is present. 
Furthermore, for objects which are small (compared to the extent 
of the PSF) the light contamination in the halo can be expected 
to be larger, in particular when the entire PN, with halo, fits on 
one IFU. In order to find out to what degree such observations 
are affected it is possible to use a procedure similar to the one 
carried out above. 



4.2. Sky subtraction 

An important issue, especially when working with weak lines 
and a relatively low resolution, is the subtraction of telluric 
(sky) emission lines from the object data. This is of special con- 
cern here since no dedicated frames of the sky were sampled, 
which have a high enough signal-to-noise in each spaxel (with 
exception of those exposures that are observed in BSW mode). 
Specifically it is the telluric Hg /14358 and the - in many circum- 
stances very - weak [O m] A4363 lines that are of main interest; 
this oxygen line is crucial when determining the electron tem- 
perature. We considered two approaches of sky subtraction. 

In the first approach an averaged sky spectrum, which is 
sampled in regions where no halo emission lines are found, can 
be subtracted from all target spectra. On the downside such a 
sky subtraction is stiff and does not account for minor differ- 
ences in line shape and wavelength calibration across the IFU. 
Consequently, for strong telluric lines residuals of "negative" 
emission can remain on one side of the line, and a residual of 
reverse sign on the other side. Such residuals in many cases pre- 
vent a correct estimate of intensities of adjacent weak lines. In 
comparison, a non-binned sky subtraction requires a separate sky 
spectrum exposure for each IFU pointing, which is observed us- 
ing the same exposure time as the object frame in order to not in- 
troduce additional noise. This is an alternative but also very time 
consuming operation. As is the nod-shuffle (BSW) technique (cf. 
Cuillandre et al. 1994; Roth et al. 2004, and references therein), 
that was used to create some of the data presented here. 

In a second approach object line fits are in every case, but 
one ([O III] /i4363), relying primarily on an appropriate fit of the 
continuum. Determining the intensity of [O iii] /14363 additional 
care is necessary, and the following procedure was adopted. 

The integrated sky spectrum of NGC 6826 (using the out- 
ermost spaxels of tile 2) is shown in Fig. 3a with an additional 
(Gaussian) line fit of Hg /i4358 in the wavelength range 4320- 
4400 A. Not only is the continuum fit a bit offset in the redmost 



Wavelength [A] 
4340 4360 4380 



Wavelength [A] 
4340 4360 4380 




4340 4360 4380 
Wavelength [A] 



4340 4360 4380 
Wavelength [A] 



Fig. 3. Illustration of the sky subtraction procedure described in 
Sect. 4.2, using data of NGC 6826. The four panels show the surface 
brightness (SB) in the wavelength range 4320-4400 A of: a) a sky spec- 
trum with a line fit of Hg A 4358, b) like a, but with ten additional lines 
that together form a set of 11 lines which purpose is to fit the sky emis- 
sion profile, c) strong object lines together with the previously fitted set 
of 11 sky emission lines, and d) that is like c, but in the much weaker 

halo. Lines starting from SB= 0.0 erg cm^^s"' arcsec^^A show individ- 
ual fits. In panels c and d the fit drawn in black is the line [O iii] /i4363. A 
total line fit is drawn on top of the observed data. Error bars are shown, 
but they are very small. 



region, but a wing of emission remains for wavelengths redwards 
of the mercury line. This wing is too wide to fit with any sin- 
gle emission line, and is likely the product of several very weak 
telluric lines; also compare with the adopted sky spectrum of 
Middlemass et al. (1991, hereafter MCWH91, see the middle 
panel of fig. 3 therein). 

In order to achieve a better fit of the continuum, and simul- 
taneously "remove" the emission wing, ten additional emission 
lines are added (at fixed wavelengths) to fit the sky spectrum. 
These additional lines are at first fitted together with the mercury 
line in a sky spectrum of each object; the result for NGC 6826 
is shown in Fig. 3b. Thereafter the intensities of these additional 
lines are fixed relative to the (variable) intensity of Hg/14358'. 



' The literature, including for example Sanchez et al. (2007) for Calar 
Alto, does not provide any spectra of the sky with a high enough reso- 
lution and at a sufficiently high signal-to-noise to discriminate between 
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All subsequent object line fits in this wavelength range (of Hy 
and [O m] /14363) are then made using this set of 1 1 lines, in 
addition to the two object lines. Provided the continuum fit is 
appropriate, and the sky spectrum is guaranteed not to contain 
any flux from the object. Fig. 3b shows that this method is accu- 
rate in estimating the sky contribution at the wavelength of the 
oxygen line. 

Figures 3c and d show two object emission line fit examples 
for NGC 6826. In addition to the sky emission lines Hy is seen 
atA = 4340 A, and [O m] at /I = 4363 A. The spectrum in panel c 
is sampled at r = 16" (arc e in tile 1, Fig. 5), and that in panel d 
at r = 26" (arc h in tile 2). Note that both object lines are much 
weaker in panel d, where a precise measurement of [O m] /14363 
would be impossible if the sky emission would not be given an 
appropriate consideration. 

The sky subtraction method described here can be used with 
any set of object and telluric lines and spectral resolutions, for as 
long as the two components are resolved. Whether it is necessary 
to fit additional Une features, as we did with the set of 1 1 lines, 
has to be decided upon on an individual basis. 

4.3. Binning of spaxels - increasing tlie signal-to-noise 

There are several ways to add individual spectra on the IFU. In 
regions where emission lines are strong enough it is, for each 
spaxel, possible to attain a surface brightness with an associated 
error. With such a two-dimensional signal-to-noise map it is then 
possible to bin the data further using, e.g., Voronoi tessellations 
(Cappellari & Copin 2003), if necessary. With weak fines, how- 
ever, which have to be binned before surface brightnesses can 
be extracted, that approach is not possible, unless a strong line 
is relied upon to create the binning map. Moreover, when the 
signal-to-noise changes very rapidly across the IFU surface - as 
it does with strong gradients - such a binning provides only a 
small improvement. Here we wanted to achieve both an increase 
in the signal-to-noise as well as a high spatial resolution. We 
therefore used two approaches simultaneously. 

For those parts of the objects where the respective emission 
fine is strong we created a two-dimensional map of the surface 
brightness. For the weaker lines in the outer parts of the objects 
the goal was to add up spectra which sample regions where the 
physical structure is (nearly) equal. In the PNe of this study all 
fines, but [O m] AA4959, 5007, Ha, and HyS are considered weak. 
Note that weak emission lines in the halo will likely, at some ra- 
dius, become readout-noise limited, in particular at higher spec- 
tral resolutions and when spaxels are small (5 O'.'S xO'.'S). Beyond 
such radfi these lines cannot be detected, even if spectra are 
binned. The halos of the objects of this study are, furthermore, 
all spherical to a good approximation, and our aim was there- 
fore to extract circular (concentric) arcs of spaxels at increasing 
radii around a defined location of the central star (CS). As is il- 
lustrated in Sect. 5, however, real objects are typically not sym- 
metric in the central region, where this binning consequently is 
less useful. A standard deviation of the intensity in each bin is 
calculated using the two-dimensional intensity maps, but only if 
an arc comprises more than three spaxels where an intensity can 
be measured. Such a value provides an estimate of how closely 
spaxels in a specific arc sample a region of similar physical con- 
ditions - a larger value indicates a deviation from spherical sym- 
metry in the physical structure. 



individual features, which is why the ten wavelengths are chosen by a 
best fit. We used the same ten values with all objects. 



A binning map illustrating the adopted configuration of spax- 
els for the observations of NGC 7662 is shown in Fig. 4. Note 
that the arcs are three spaxels wide in the outermost part of the 
halo (to the right), in order to compensate for the lower line in- 
tensities expected there (cf. Sect. 5.1.1). Showing the diff'erence 
in the binning setup between the LARR and PPAK observations, 
the two biiming maps used in the analysis of the NGC 6826 data 
are shown in Fig. 5. All arcs in both tiles are one spaxel wide, as 
a two-spaxel wide arc was not found to provide any significantly 
better signal-to noise in this case. 

4.4. Dereddening 

We measured extinction coefficients using the Hy/H/? line ratio, 
but we also used values from the fiterature. In several cases we 
could calculate a value across the PN surface, but in the faint 
halos we always used a constant. In Sect. 5 we compare lit- 
erature and measured values for each object separately. In the 
dereddening procedure we used the interstellar extinction curve 
of Fluks et al. (1994). When measuring Hfi we did not consider 
He II emission, since it is negligible in most cases. 



5. Results 

We made all line fits using the IFU analysis package, iFSFn, that 
we developed for this purpose, ifsfit is based on the idl rou- 
tine MPFir of the freely available Markwardt library^ and was 
designed to work with most IFUs. Gaussian curves were used to 
fit emission fines, and polynomials to fit continua. 

In a first step we made a line fit for each line and spaxel on 
every data set in order to get a surface brightness and an associ- 
ated error. The outcome form two-dimensional spatial maps of 
surface brightnesses. Not all fine fits are good, and bad fits and 
inaccurate fits are marked as such after a visual inspection of 
each individual line fit. A bad fit typically occurs when a line 
disappears in noise, or if an unmasked line was removed as a 
cosmic ray by accident during the data reduction (see Sect. 3). 
Inaccurate fits may occur if the continuum level is badly esti- 
mated, or if the actual line width in a particular spaxel is larger 
or smaller than the pre-defined and allowed deviation of ±15%, 
from the theoretically determined full width at half maximum 
(of the grating in the used wavelength range). 

In a second step we used a pre-defined map of radially 
binned spaxels (see Sect. 4.3), where we averaged the contribut- 
ing spectra before fitting the emission lines. The details of the 
fitting procedure are the same as in the first step. The outcome 
consists of a surface brightness and an error for each (radial) 
position of the respective binned arc and emission line. Surface 
brightness errors are typically small. When we draw error bars in 
plots of the radial surface brightnesses we instead use the spaxel- 
to-spaxel standard deviation (Sect. 4.3). 

Concerning physical quantities an electron temperature (Tg) 
was in every case calculated using the expression of Osterbrock 
& Ferland (2006, see eq. 5.4 therein) for [Oiii], assuming a neg- 
figible dependence on the electron density (A^e)- An electron 
density was derived using two methods. Where available the 
forbidden line ratios [ Ar iv] ^^47 1 1 /4740, [On] AA3126/3129, 
and [Sii] /1/167 17/6730 were used to calculate A^e.Anv, A^e.On, 
and A^e.Sii - using the iraf nebular package of Shaw & Dufour 
(1995). We calculated maximum error limits of the density using 
the minimum and maximum estimates of the respective ratio. 



^ http://cow.physics.wisc.edu/~craignVidl/fitting.html. 
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Fig. 4. Binning map used in tlie analysis of our observations of NGC 7662. Each spaxel (square) represents an area of I'.'O x I'.'O on the sky, and 
data is here present out to a distance of 36" from the CS (•). Each of the three tiles, in the east-west direction, has 16x16 elements and represents 
the shape of the PMAS LARR IFU; in order to emphasize that each tile is observed separately they are slightly separated from each other in 
the figure. Spectra of spaxels marked with the same letter are summed up to increase signal-to-noise (cf. Sect. 4.3) - forming concentric arcs of 
spaxels, which are centered on the CS. Gray shades of the spaxels were chosen to enhance the shape of the arcs. Spaxels marked with the symbol 
® are not used. The N-E most sides of spaxels in tiles 1 & 2, together with the S-E most sides of tile 3 were all masked to not be used since they, 
due to the DAR correction procedure, depend on flux outside the observed area (see Sect. 3). Spaxels marked with the sjmbol ♦ were not used 
since either one or several emission lines were removed as cosmic rays in the data reduction; note that the DAR correction procedure causes each 
non- wanted spaxel to result in four masked spaxels. The ♦ symbol is also used to mask most of tile 1, where the object structure is asymmetric. 
This map can be compared with the object image in Fig. la, and also the line intensity maps in Fig. 7. 
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Fig. 5. Biiming map used in the analysis of the observations of NGC 6826. Note that north is down and east to the right. In difference to the LARR 
IFU, which is arranged in a square structure using square-shaped lenslets, the PR\K IFU is arranged in a hexagonal structure where each fiber is 
circular. Each spaxel (circle) represents an area of 5'.' 64^ on the sky, and the outer tile (2; right) reaches a radius of 80" (in arc w). The coordinates 
indicate where the center of the respective tile (marked with a white upside down triangle) is with respect to the central star (that is indicated in 
the left tile with a bullet). These maps can be compared with the object image in Fig. le, if they are flipped in both directions. For further details 
of the figure see the caption of Fig. 4. 



Since the (electron) density in the halo in general is so low 
that forbidden emission line ratios always are outside density 
sensitive ranges (SK89) another method must be used there. 
Assuming spherical symmetry, full ionization, and a distance it 
is possible to derive an emissivity (j) using an (inverse) Abel 
transform on the radial surface brightness structure of a recom- 
bination line (such as Hy6). The electron density is proportional 
to the emissivity as Ne^c f^^. This method is first used by Plait 
& Soker (1990, hereafter PS90) on NGC 6826 in order to calcu- 



late an electron density structure. The electron density can there- 
after be converted into a mass density and a mass loss rate (see 
Sect. 6.2). Since our density estimates are based on the observed 
intensity structure, and not on the true (deconvolved) intensity, 
they are partly affected by scattered light, which is the strongest 
in the inner halo (cf. Sect. 4.1). 

We present the results for each object separately. In consec- 
utive order they are surface brightnesses, electron temperatures, 
and electron densities. 
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Fig. 6. A map of the calculated extinction coefficient c of NGC 7662 for 
tiles 1 and 2. The CS position is marked with a filled star; the two num- 
bered circles indicate regions observed by Barker (1986). On the east 
and north sides the "missing" column and row of values were removed 
in the DAR correction procedure. The vertical dash-dotted line marks 
the tile boundary. For further details see Sect. 5.1. 



5.1. NGC7662-PNG106.5-17.6 

The spaxel binning map that we used with the data of this ob- 
ject to calculate a radial intensity structure is shown in Fig. 4. 
In the central region NGC 7662 is asymmetric, which is why 
most of the spaxels there were masked, in order to achieve a 
smoothly varying surface brightness across the inner-middle tile 
(1-2) boundary. A recent detailed discussion of structural prop- 
erties of NGC 7662 is given by Guerrero et al. (2004). Moreover, 
a distance is determined by Zhang (1995, d = l.ll-1.24kpc, 
who includes a literature summary), Hajian & Terzian (1996, 
d = 0.79+0.75 kpc), and CSSP03 {d = 2.2 kpc, from Gorny priv. 
comm.). MCWH91 use d - 1.5 kpc. We used d = 1.7 kpc (see 
Sect. 6.2). A value of the effective temperature of the CS is given 
by Gorny (priv. comm.), retf(Zanstra) = 100 000 K. 

Numerous references give one (or a few) scalar extinction 
coefficient for this extended object, see, e.g., Osterbrock (1964, 
c = 0.42), Peimbert & ToiTes-Peimbert (1971, c = 0.16 + 0.12), 
Harrington et al. (1982, c = 0.23), Kaler (1986, c = 0.17, hereafter 
K86), MCWH91 (c = 0.23), Lame & Pogge (1996, c = 0.16 ± 
0.04), and Hyung & Aller (1997, c = 0.10). We first calculated 
a map for the central region of the nebula, see Fig. 6. Two out 
of six positions where Barker (1986, qi) = 0.26, C(2) = 0.17) 
measures the extinction coefficient are also shown in the map, 
indicating a fair agreement. A vague outline of the rim is also 
visible, where the extinction is slightly higher than elsewhere. 
Comparing this map with that presented in Sandin et al. (2007, 
fig. 3a, that uses the same data, but does not correct for DAR), it 
is seen that there is less structure in this case. The average of the 
values in the central tile is c = 0.24 + 0.09, which is close to the 
value measured by MCWH91. In the outer parts of the nebula, 
where an extinction coefficient cannot be determined (as in the 
westmost side of Fig. 6), we always use this constant value. 

5.1.1. Surface brightnesses 

Two-dimensional surface brightness maps of four selected emis- 
sion lines are shown in Fig. 7. The (asymmetric) elliptical rim 
is clearly seen (in the central PN, i.e. the leftmost tile) along 
with a rapidly radially outwards decreasing intensity. Each panel 
shows contours of the respective logarithmic intensity, revealing 
the geometrical shape of the nebula at different intensity lev- 
els. In the two outer tiles intensity contours (which are circular 



to first order) are present almost all the way to the west-most 
(right) boundary in the strong oxygen line [O m] /15007 (Fig. 7a). 
H/3 (Fig. 7b) shows larger fluctuations in the halo already in 
the middle tile. Also note that this recombination line is struc- 
tured differently on the rim compared to the other three lines; 
the rim is closed on the southeast side of the CS. [Oiii] /14363 
and [Ne m] /i3869 (Figs. 7c and d), finally, are both weak and 
were not detected either in the outer part of the middle tile or in 
the outermost tile. 

Corresponding signal-to-noise maps are shown in Fig. 24 
for the two oxygen lines. The discontinuous contours, that are 
seen at tile boundaries in Fig. 24, indicate a changing signal- 
to-noise with exposure time. Such a discontinuity is not visible 
between the outer two tiles of the [O m] A4363 line that is shown 
in Fig. 24b. In this case no line could be measured in the outer 
parts of the middle tile. Additionally, note the intensity fluctu- 
ations at lower signal-to-noise values (S/N ^ 10) as seen in the 
rightmost tile (3) of Fig. 7b and Fig. 24b. Similar intensity maps 
were created for all lines and objects, but are not all presented 
due to limitations of space. 

In Fig. 8 we present radial surface brightness structures that 
were derived from the binned arcs for 13 emission lines. Note 
that we also could measure the unresolved oxygen line doublet 
[On] AA3126, 3729 in the inner region of the central tile. Most 
lines drop in intensity by a factor of 3000 and more between 
the central parts and the outermost parts of the observed halo. 
The argon and helium lines are not present in the region around 
r w 21", since the exposure time used with tile 2 (200 s) did not 
allow a detection in these lines. Due to the asymmetric shape of 
the central nebula, line intensities in adjacent spaxels of tiles 1 
and 2 differ by as much as 50 per cent, explaining the rather large 
step in intensities at r ^ 8". Likewise, signs of asymmetry are 
also visible for intensities at r ^ 24" (arc f) in tiles 2 and 3; this 
arc comprises five and ten spaxels, respectively. The size of the 
error bars in the right hand side panels show that the departure 
from spherical symmetry in separate bins is the largest in all lines 
for6"$rg8" and 12"<r$16". 

Comparing the PN proper and the halo further the relative 
intensities (right hand side panels in Fig. 8) show an increas- 
ing trend in the halo outwards from the shell-halo transition, i.e. 
for r ^ 13", in all (but the hydrogen) lines. The intensity ratio 
Hy/HyS (HS/Hf)) is, where a two-dimensional extinction coeffi- 
cient map could be used, constant to a good approximation. In 
Case B ionization the, relatively temperature insensitive, ratio is 
0.47 (0.26) for = 10 000 K (see, e.g., table 4.4 in Osterbrock 
& Ferland 2006). While the Hy/HyS-ratio is close to the theo- 
retical value the H(5/H/3-ratio is off by about 15 per cent out 
to the shell-halo transition. This difference indicates a system- 
atic error in the surface brightness of H6, the agreement between 
theory and observations of this ratio is better with the other ob- 
jects. Emission lines are in every case, but H6 and He 1/14471, 
present at intensities below 1 x 10^'^ergcm^^s^'arcsec"-. The 
sky continuum intensity is about 1 x 10"'^ ergcm^^s^'arcsec"^, 
which is the sensitivity limit achieved with the adopted instru- 
mental setup and exposure times^. The amount of scattered light 
is small, this is confirmed when comparing surface brightnesses 
with the radial intensity level in fig. 21 of CSSP03 ([Oiii] and 
[O III] - diffuse light, where it was necessary to subtract scattered 
light present in JKT data). 



^ For the bluest measured emission line, [Ne iii] /i3869, the sky con- 
tinuum intensity is about 2x10"'^ ergcm^^s^'arcsec"^. 
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Fig. 7. Surface brightness maps of a mosaic 
of the three observed tiles of NGC7662 for 
four emission lines. Intensities are logarithmic 
and the unit is ergcm"-s"'arcsec"-. The verti- 
cal dash-dotted lines mark the separation be- 
tween the three tiles. The position of the cen- 
tral star is marked with the symbol Colors 
repeatedly show additional contour lines; the 
separation between two color contours corre- 
sponds to a 0.083 dex = 21% change in sur- 
face brightness. Dense lines in the northmost 
row of spaxels, in the middle tile, are cre- 
ated by the contour drawing routine, and do 
not show a real gradient. From the top the 
panels show four increasingly blue lines: a) 
[Oiii]/15007, b) HyS, c) [Oiii]/l4363, and d) 
[Neiii]/i3869. Corresponding signal-to-noise 
maps for [Oiii] .15007 and [Ora]A4363 are 
shown in Fig. 24. The dotted line in the mid- 
dle tile indicates the location of the shell-halo 
transition (rsheii = 13"). For further details see 
Sect. 5.1.1. 



5.1.2. Electron temperature 

We could measure [O iii] /14363 out to a distance of r = 30" from 
the CS. In the calculation of the temperature we used all three 
diagnostic [O iii] lines. The radial structure is presented in Fig. 9. 
Additionally, we draw two gray lines in the figure, that were 
derived using a constant extinction coefficient; c = 0.1 (0.4) was 
used for the lower (upper) line. These two lines demonstrate how 
sensitive the measured temperature is to the use of the extinction 
coefficient. The difference in the halo (for r ^ 13"), where the 
value of the actual extinction coefficient is more uncertain, is 
about 500 K between the black line and the lower gray line. 

From the shell-halo transition, and outwards in the halo, the 
temperature increases from = 12 500 + 100 K, at r = 13", to 
Te = 18 500 + 1000 K, for r > 22". Moving inwards, instead, the 
temperature also increases, from Tg = 12 500 K, at r = 13", to 
about Te - 14 700 ± 300 K at the location of the rim; due to the 
asymmetric inner region this increase is not smooth (see below). 

The halo measurement of MCWH91, Te = 17 500 + 1700 K, 
is also indicated in Fig. 9. It was made in about the same region 
of the halo that is covered here. The agreement with our value at 



r = 30" is fairly good. Additionally, MCWH91 quote an average 
core temperature of = 13 1 00+ 500 K (where the measurement 
of [O III] /14363 is taken from Peimbert & Torres-Peimbert 1971). 
Hyung & Aller (1997) derive Te ^ 12 500K (using c = 0.10 and 
assuming A^e - 8000 cm"-') in about the same region as position 2 
of Barker (1986). Further estimates are presented by, e.g., K86, 
Peimbert et al. (1995) and Zhang et al. (2004). The values are 
in general about lOOOK lower than what we found here in the 
central nebula. 

In comparison to the binned radial temperature structure we 
show a two-dimensional temperature map of the central regions 
in Fig. 10. A comparison is best made with the central parts of 
the map shown in fig. 4 of Reay & Worswick (1982, that is 
derived using an electronographic method). The agreement is 
good - the asymmetric region extending to the west is clearly 
seen in both maps. Absolute values are also in better agreement 
than with the other references. Barker (1986), moreover, derives 
a temperature at five different positions across the central neb- 
ula, where two are covered in Fig. 10 (positions 1 [4"N] and 2 
[3'.'7 E, 5" S]). Both values, 13 800 + 500K and 13 100 + 500 K, 
respectively, are again about 1000 K lower than what is shown 
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Fig. 8. Radial surface brightness structures of NGC 7662, cf. Sect. 5.1. All ordinates are logarithmic. Intensities of each ion are indicated using the 
same symbol in both left and right hand side panels. Panels on the right hand side show the log of each intensity ratio relative to H/J. From the top 
the panels show 13 emission line intensities of five elements: a) hydrogen, b) helium, c) oxygen, and d) neon and argon. Error bars due to Poisson 
and readout noise were omitted since they in every case are smaller than the symbol sizes. Instead error bars indicate the spaxel-to-spaxel standard 
deviation - for clarity only for a few lines (cf. Sect. 5). The radial positions of the three tiles are indicated in the upper part of panel d. For each tile 
the inner edge is marked by an upwards pointing marker and the outer edge by a downwards pointing marker. Note the overlap in the data between 
tiles 1 & 2 and 2 & 3 (see Fig. la for the actual positions of the tiles on the sky). Vertical gray bars indicate the location of the rim at rnm = 6", and 
the shell-halo transition at rjheu = 13". The horizontal dashed gray lines in the right a panel indicate line ratios for Case B recombination theory at 
Te a 10000 K. Horizontal dash-triple-dotted gray lines in the right hand side panels are guides, that indicate the constant ratios 0.01, 0.1, and 1.0. 
Data is missing for some of the weaker lines at r a;20", since the exposure time of tile 2 (200 s) was not long enough to measure them accurately. 
The offset between data of tiles 1 and 2 at r « 9" (and tiles 2 and 3 at r » 24"), is due to deviations from symmetry. 



in Fig. 10. In the halo the low signal-to-noise in each spaxel of 
the weak [O m] /14363 prevents the calculation of a map, instead 
it is necessary to bin the spaxels, as we did for Fig. 9. 

5.1.3. Electron density 

At first we measured the electron density at a few radial loca- 
tions using the [ Ar iv] /147 1 1 /4740 and [O n] ^3726/3729 line 
ratios. Using the oxygen doublet we derived ratios for r = 



8'.'5, 9'.'5, and \Q(!5 to find N^^o^ = 1800:1:^^^, 1540;!:^^, and 
1260^[jQ cm"^. Due to the low resolution it is difficult to deblend 
[Ariv]/I4711 and the, in this case, very weak He 1/14713. The 
resulting values for r = 8'.'5, 9'.'5, KX'S, ll'.'S, and 12'.'5 were, 
A^e,Ariv=4540 ± 230, 4110 ± 200, 3950 ± 170, 3760 ± 160, and 
3460 ±180 cm~^. Our values are in good agreement with Lame 
& Pogge (1996), who measure A^e.Sn = 4000 cm"^ in the shell. 
Other authors report on values measured in regions not consid- 
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Fig. 10. Electron temperature map of the cen- 
tral regions of NGC7662, shown in the unit 
of 10^^ K. Colors repeatedly indicate additional 
contours at a 200 K interval. The vertical dash- 
dotted line marks the separation between tiles 
1 and 2. The position of the CS is marked with 
a filled star Temperatures are measured in the 
interval 1 1 800-16 000 K. The location of the 
shell-halo transition is marked with a dotted 
arc, Tshcii = 13". The lowest values are found in 
the outer shell and the highest values in regions 
inside the rim, compare with the surface bright- 
ness maps in Figs. 7a & c. The two numbered 
circles indicate regions also observed by Barker 
(1986). For further details see Sect. 5.1.2. 
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Fig. 9. Radial temperature structure of NGC 7662. In addition the halo 
measurement of MCWH91 is indicated with a triangle (A). The inner 
gray vertical line indicates the location of the rim, and the outer line the 
shell-halo transition. The upper (lower) gray line shows calculated 
with the extinction coefficient set to c = 0.4 (0.1). This radial structure 
can be compared with the two-dimensional temperature map of the in- 
ner parts of the nebula in Fig. 10. For further details see Sect. 5.1.2. 



ered here, e.g.: Barker (1986, A^eAr,v(l,2) = 6300, 7000 cm"^ 
and A^e,ciiii(l , 2) = 3400, 1500 cm-^. 

On the other hand we calculated a density structure using the 
inverse Abel transform with H/3. The resulting electron density 
structure is shown in Fig. 11, together with the four measured 
values of the argon line ratio, against which we scaled the density 
structure. In drawing the upper axis we used the distance d - 
1.7 kpc. Note that the density is less accurate for the central and 
outer regions (10" ^ r ^ 30") due to the asymmetric center and 
an H/3 structure, which is truncated at r = 35" (since we did not 
cover the full radial extent of the halo). 
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Fig. 11. Radial density structure of NGC 7662 (diamonds, 0) calculated 
using an Abel transform. The structure is shown with a logarithmic or- 
dinate in panel b (bottom). In panel a (top) our mass loss rate estimate 
M/u is also shown for a part of the radial extent (using squares and 
the two right axes, cf. Sect. 6.2). The gray axis shows the mass loss 
rate normalized with the maximum plotted value (M„,ax), in per cent. 
The temporal axis above the plot shows the age of the structure radi- 
ally, assuming a constant outflow speed of m = 1 km s" ' . The thick gray 
line shows a fit of the density to the power law A'e oc r^" for the region 
17" < r < 22". Bullets mark individual measurements of the electron 
density A'c.ahv, along with maximum error estimates. The inner gray 
vertical line indicates the location of the rim, and the outer line the 
shell-halo transition. For further details see Sect. 5.1.3. 



5.2. IC 3568 - PNG1 23.6+34.5 

IC3568 is a close to round nebula with a measured size of 
'"shell = 9" (Collins et al. 1961, who used observational data of 
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Capriotti & Daub 1960), or rsheii =6" (see Zhang & Kwok 1993, 
and references therein). Distance determinations, including lit- 
erature summaries, are made by, e.g., Zhang (1995, d - 2.68- 
3.12kpc) and Phillips (2005, £/= 1.08-3.91 kpc). CSSP03 quotes 
a Shklovsky distance of d - 4.5 kpc, we used d - 4.3 kpc 
(see Sect. 6.2). The effective temperature of the CS is deter- 
mined to be Tetf = 50 000 K (Harrington & Feibelman 1983), 
Teff = 51 300K (Preite-Martinez et al. 1991), and reff(Zanstra) = 
55 000 K (Gorny priv. comm.). 

Values of the extinction coefficient are given by, e.g., Collins 
et al. (1961, c = 0.25), Peimbert & Torres-Peimbert (1971, c = 
0.10 + 0.12), Sistla & Kaftan-Kassim (1976, c = 0.23+0.05), 
Pottasch et al. (1977, c = 0.29+0.04), Barker (1978, c = 0.252 + 
0.032), Harrington & Feibelman (1983, c = 0.25), K86 (c = 0.19), 
Tylenda et al. (1992, c = 0.4), Kwitter & Henry (1998, c = 0.45, 
measured 4" N of the CS, hereafter KH98) and Liu et al. (2004, 
c - 0.26), who all report on values which are measured at or 
close to the nebular center. Our measurements in the central tile 
revealed a varying coefficient, that is c ^ 0.4 4" west of the CS, 
c^O.2 4" north of the CS, and negative values are found east of 
the CS. Instead of using a map we chose to set c = 0.2 across the 
entire surface of the nebula. Also Phillips et al. (1979) present 
evidence for a differential extinction. 

The two binning maps that we used with the data of this 
object are illustrated in Fig. 25. These maps are identical with 
the exception of the DAR mask, which is different for wave- 
lengths blue- and redwards of the selected reference wavelength 
(/lief = 5100 A; cf. Sect. 3). Data are here present out to a distance 
of r = 25" from the CS. 

5.2.1. Surface brightnesses 

Radial surface brightness structures, which were derived using 
the binned arcs, are presented in Fig. 12. From the spatial struc- 
tures of the lines the location of the rim and the shell-halo tran- 
sition were determined to be riim - 2'.'5 and rsheii - 9'.'0, respec- 
tively. Physical properties differ in the spaxels of those arcs, that 
are present in both tiles 1 and 2, causing a step of :$ 0.2 dex in 
the surface brightness structure for 9'.'0 ^ r 12". In addition 
to the presented lines [Nii] /16548 and [Ariv] AAIW were also 
detected, but could not be accurately determined. The nitrogen 
line is very weak compared to Hot, and the argon line cannot be 
deblended from He i /14713 due to the low resolution. 

Hy6, the two oxygen lines [Om] /1/14959, 5007 and neon 
[Neiii] /13868, are all present in the entire region covered. The 
surface brightnesses in the outermost arc in the halo are about 
0.6-2 X 10"^ times as strong as on the rim. Also note that the 
shell, compared to the rim, is weaker than for the other objects; 
the corresponding intensity ratios are about 25, instead of about 
5-8. Two out of three of the relevant halo selection criteria of 
CSSP03 are met. The first one is not, since we did not cover the 
full radial extent of the assumed halo. Since the [O m] emission 
is too strong in the outer regions, to be caused by scattered light, 
we conclude, in contrast to the result of CSSP03, that IC 3568 
has a halo. 

5.2.2. Electron temperature 

We calculated a temperature using all three diagnostic [Om] 
lines. A two-dimensional map of the central region is shown in 
Fig. 13, revealing a slightly asymmetric structure. We measured 
temperatures in the interval 10400 % % 12 000K. A simi- 
lar temperature map is presented by Reay & Worswick (1982, 
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Fig. 13. Electron temperature map of the central regions of IC 3568, in 
the unit of 10^ K, cf. Sect. 5.2.2. The location of the rim is marked with 
the inner dotted circle, and the shell-halo transition with the outer circle. 
For further details of the figure see Fig. 10. 

10 800 < Te < 13 200 K). A comparison between the two maps 
shows few similarities, although there is a match in the existence 
of a region with lower temperatures east of the CS. 

We also derived a temperature structure from the radially 
binned data, that is shown in Fig. 14, including the inner halo; 
error bars are notably larger in the halo than in the central neb- 
ula. The choice of concentric arcs to bin the spaxels in the region 
inside the shell is, as the two-dimensional map in Fig. 13 shows, 
not optimal (but better than for NGC7662). Compared to the 
temperature structures of NGC 7662 and NGC 6826 (see Figs. 9 
and 20) there is a weaker correlation between the temperature 
minimum at r ^ 6" and the location of the shell-halo transition 
('"shell - 9"). IC 3568 also, however, shows an increasing temper- 
ature gradient in the region 6" < r < 12", where the temperature 
increases from Te = 1 1 000 K to ^ 14 000 K. 

Previous derivations of the temperature in the central nebula 
are presented by Barker (1978, = 10 800+500 K; with revised 
values calculated by K86, - 10 410 K, and McKenna et al. 
1996, Te = 9990 K), KH98 (reo„, = 10500 + 5000 K, TeNu = 
6900 + 5000 K), Zhang et al. (2004, Te^o,,, = 11 370 + 200 K and 
T^e.Ba = 9300+900 K). Most of these values are lower than what 
we calculated for the central nebula. The temperature map of 
Reay & Worswick meanwhile shows a larger value in the center 
(by about 1000 K). 

5.2.3. Electron density 

A review of electron density measurements in the literature 
shows large differences: Barker (1978, A^eOn - 10000 ± 
9000cm-3), SK89 (A^e,o„ = 3390cm-3 and A/^e,s„ = 3090cm-^ 
that were derived using observations of previous authors), KH98 
(measure A^e s i. = 900+700 cm^^ 4" north of the CS), Wang et al. 
(2004, A^e,o„ = 2UQ^%cm-\ N,,su = 1550!«^Ocm-3, N,,m,, = 
1510!23°cm-^ and A^e.cim = 39+22" cm^^), ^j^^ 2hang et al. 
(2004, logA^e.Bai = 3.8+0.2). Our observations of [Ariv] ^4711 
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Fig. 12. Radial surface brightness structures of 1C3568, cf. Sect. 5.2.1. From the top the panels show 16 emission lines of seven elements: a) 
hydrogen, b) helium, c) oxygen and nitrogen, and d) neon, argon, and sulfur. For further details of figure properties see Fig. 8. Vertical gray bars 
indicate the location of the rim at ^ 2'.'5, and the shell-halo transition at rsheii - 9". Emission lines with wavelengths longer than Ata = 5 100 A 
are drawn with gray lines, remaining lines are drawn in black. Also see Sect. 5.2. 



are too inaccurate to make a density estimate - the actual density 
also appears to be below the conservative density sensitive range 
of [Ariv] (SK89; logATe.Anv^S.S). 

Instead we calculated a radial density structure using H/3, 
with the assumptions of A^e(core) ^ 2000 cm"-' and d = 4.3 kpc, 
see Fig. 15. This density is less accurate for the outer regions 
(r ^ 20"), since we did not cover the full radial extent of the 
halo. 

5.3. M2-2 - PNG147.8+04. 1 

Existing studies of this small and distant nebula mostly discuss 
its geometrical shape and distance. For instance, Cahn et al. 
(1992, and references therein) measure a diameter of the inner 
nebula of T.'O. Zhang (1995) finds the same size and derives 
(and quotes) distances of d = 3.9^.4 kpc. Moreover, GVM98 



(Manchado et al. 1996) find the size of the inner nebula to 
be 5'.' 3 X 4'.'0, and measure a halo in Ha of the proportions 
17'.'4 X 14'.'9. An [O iii] image of the nebula is provided by Chu 
(1989). AUer & Keyes (1987), finaUy, calculate an elfective tem- 
perature Teff = 80 000 K of the CS. 

Values of the extinction coefficients are derived by Kaler 
(1983a, c = 1.25 ±0.15), Kaler (1985, c = 1.58), Aller & Keyes 
(1987, c = 2.0), and Tylenda et al. (1992, c = 1.46). Using the 
same procedure as with NGC 7662 we first calculated a map of 
the extinction coefficient in the central tile. The average of this 
map is c = 1.50 ± 0.25, that agrees well with the literature val- 
ues. The map is used to deredden the data in the blue wavelength 
range for the central tile, before binning. In comparison we used 
the average value in tiles 2 and 3, and for aU data in the red 
wavelength range. 
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Fig. 14. Radial temperature structure of IC 3568, cf. Sect. 5.2.2. For fur- 
ther details of the figure see the caption of Fig. 9. 
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Fig. 15. Radial density structure of IC 3568, cf. Sect. 5.2.3. The thick 
gray line shows a fit of the density to the power law A'e oc r"" for the re- 
gion 10'.'5 < r< 12'.'5. For further details of figure properties see Fig. 11. 

The two spaxel binning maps that we used in the blue and 
red wavelength ranges are illustrated in Fig. 26. These maps 
differ due to the reverse shift of the DAR masks for wave- 
lengths shorter and longer than the reference wavelength /lief - 
5050 A. Several emission lines, such as [Nii] AA654S, 6583 and 
[S ii] /1/16717, 6730, are predominantly present in ansae on the 
south side of the nebula (see the inset in Fig. Ic and also fig. Id 
of GVM98). In order to not affect the halo measurements when 
binning the spaxels we masked these regions in tiles 1 and 2. 



5.3.1. Surface brightnesses 

Radial surface brightness structures, that we derived using the 
binned arcs, are presented in Fig. 16. We determined the lo- 
cations of the rim and the shell-halo transition from the radial 
structures of the hydrogen and [Om] lines, as Trim - 2'.'0 and 
'"shell - 6'.'5. The two emission lines [Om] AA4959, 5007 were 
present out to a radius of r = 17" and r = 20" from the CS, respec- 
tively. In comparison, H/3 was detected out to r = 8'.'5, beyond 
this radius this line is too weak to be detected. The resulting halo 
diameter of D = 40" is consequently 2.3-2.7 times larger than the 
halo size reported by GVM98 (as specified above) and also 2.7 
times larger than the diameter Phillips & Ramos-Larios (2006) 
measure in the J band (D = 15", see fig. 6 therein). Furthermore, 
all intensity ratios we measured here between r = 8"0 and the 
bright center are 25-100, which are a factor of ten larger (and 
more) than the ratios of 2-5 found by GVM98. Finally, from a 
scattered light test, that we made following the method described 
in Sect. 4. 1, and using [O iii] /15007, we find graphs that are very 
similar to those shown in Fig. 2 for NGC7662 (although with 
the smaller radii of M 2-2). The measured halo of M 2-2 is thus 
real. 

Showing the, to a good approximation, spherical symmetry 
of the structure two-dimensional surface brightness maps of four 
selected emission lines are illustrated in Fig. 27. The ansae, men- 
tioned in Sect. 5.3, are not clearly visible in the presented lines, 
with the exception of the slightly increased value seen in the S-E 
corner of the middle tile (right hand side panels in the figure). 

5.3.2. Electron temperature 

We calculated a temperature using all three diagnostic [Oiii] 
lines, where the auroral line [O iii] /14363 could be measured out 
to a distance r = 7"25, i.e. barely reaching into the halo. The 
resulting radial structure is presented in Fig. 17. Note that the 
temperature is very high in the center, ^ 20 000 + 900 K, and 
much lower (and more or less constant) in the remaining parts, 
000 + 200 K. This average value is in perfect agreement 
with the value of Aller & Keyes (1987, = 12 000 K). Using 
infrared observations Phillips & Ramos-Larios (2006) note that 
the index ( J - H)q decreases systematically away from the CS 
of M2-2. The authors discuss the possibility to explain such a 
decrease in terms of hot dust continua. There is in the presented 
data no clear indication of an increasing temperature in the halo. 

5.3.3. Electron density 

Due to a comparatively low spectral resolution of the obser- 
vations we could not accurately deblend the adjacent lines 
[ Ar iv] AAl 1 1 and He i AAl 1 3 (where the helium line is found 
to be the weaker). Instead we measured the sulfur line ratio 
[S ii] AA6111 16131 at r = 5'.'75, 6'.'25, and 6'.'75, and found the 
ratios 0.698 + 0.059, 0.728 + 0.066, and 0.829 + 0.073, respec- 
tively. These ratios give the densities A^e.Sn =2420+^2„, 2080+^^°, 
and 1310^34Qcm"-^. In comparison Aller & Keyes (1987) mea- 
sure A'e.sn = 5000 cm"^ in the central nebula. Additionally we 
calculated a radial density structure using the Abel transform 
with Wa, that we measured throughout most of the halo. When 
scaling the density structure we assumed a distance of (i = 4kpc. 
As for NGC 7662 and IC 3568 we did not measure the outermost 
part of the halo, which is why the density is more uncertain for 
r ^ 16". Due to asymmetry in the Ha-structure we did not calcu- 
late a density in the central region either. Both the radial structure 
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Fig. 16. Radial surface briglitness structures of M2-2, cf. Sect. 5.3.1. From tlie top the panels show 18 emission lines of seven elements: a) 
hydrogen, b) helium, c) oxygen, nitrogen, and sulfur, and d) neon and argon. Data of lines drawn in black in the blue wavelength range were 
observed 04. 1 0.2007, and data of lines drawn in gray were observed 05. 10.2007. Vertical gray bars indicate the location of the rim at ^ I'.'O, and 
the shell-halo transition at Tsheii - 6'.' 5. The horizontal dashed gray lines in the right hand side panel a indicate line ratios for Case B recombination 
theory at « 12 000 K. For further details of figure properties see Fig. 8. 



and the individual measurements are presented in Fig. 18, along 
with the mass loss rate (cf. Sect. 6.2). 

5.4. NGC 3587 -PNG1 48.4+57.0 

NGC 3587 has a thin halo outside an extended bright region. The 
diameter of the PN proper is D > 200", and surface brightnesses 
are low. Moreover, since the effective temperature of the CS is so 
high(reff = 93 900+5600 K, Napiwotzki 1999; Teff^ 110000 K, 
McCarthy et al. 1997), and the luminosity so low (L = 41- 
148 Lq, Kaler 1983b; L= 63 L© Gorny priv. comm.), NGC 3587 
is suggested to have a so-caUed recombination halo (Tylenda 
1986; Corradi et al. 2000; CSSP03). The halo, that has an asym- 
metric shape, suggests an interaction with the surrounding inter- 
stellar medium (see, for example, Sabbadin et al. 1985; Guerrero 



et al. 2003, and references therein; the latter reference also pro- 
vides a list of distance measurements). We observed the shell- 
halo transition region east of the CS in order to find out what 
the detectable radial extents of various high- and low-excitation 
emission lines are, and also to measure a temperature. The spaxel 
binning map we used with the observations in the red wavelength 
range is shown in Fig. 28. We did not calculate an extinction co- 
efficient, but instead used the constant value c = 0.0 (see, e.g., 
Kaler et al. 1990; Tylenda et al. 1992). 

5.4.1. Surface brightnesses 

We present radial surface brightness structures, which were cal- 
culated using the binned arcs, in Fig. 29 for 16 observed emis- 
sion lines. Observational conditions allowed lines in the red 
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Fig. 18. Radial density structure of M 2-2, cf. Sect. 5.3.3. The thick gray 
line shows a fit of the density to the power law A'e r^", for the region 
7'.' 5 < r < 9'.'5. Bullets with error bars show individual measurements of 
^cSii- For further details of figure properties see Fig. 11. 



wavelength range to be measured in a larger spatial domain, 
reaching inwards to r = 65". The agreement between data from 
the two different years is good, small deviations should be due 
to sampling in slightly different regions (see Fig. Id). Compared 
to the other objects of this study the intensities decrease across 
the observed region by less than a factor of 100. The shell-halo 
transition at r ^ 105" is less sharp than for other objects (also 
compare with the east [left] side of the two images in fig. 10 of 



CSSP03). The Hne intensity structures of Ha and [Oiii] /15007 
agree well with fig. 19 in CSSP03 - the rim and shell-halo tran- 
sition locations were taken from this reference. Furthermore, as 
is seen in Fig. 29a (right hand side panel) the agreement with 
recombination theory is fair for Ha and Hy where r < 105". H5 
is more noisy. Uncertainties are much larger for all lines in the 
halo. 



5.4.2. Electron temperature 

K86 and McKenna et al. (1996) measure an [Oiii] electron tem- 
perature of Te = 10910K and Te = 10 840K, respectively. Both 
papers apply updated atomic data on the observed line ratio of 
Torres-Peimbert & Peimbert (1977, who put the slit on top of 
the CS). Kwitter & Henry (2001) measure TeOm = 10 600 K, 
re,N„ = 9400 K, re,o„ = 11 600 K, and Te.sm = 12400K (all 
+ 10%). Using the values of [Oin]/l4363, in the radial range 
85" < r < 95", we found a more or less constant temperature, 
which average value is re^om = 11 500 + 240 K. A comparison 
with the values of the central region indicates a close to constant 
temperature across the nebular surface. 



5.4.3. Electron density 

Two of the line doublets we observed are indicators of the elec- 
tron density, viz. [On] and [Sii]. However, the density in the 
shell and halo is low. So low that the measured line ratios of 
these two indicators are below the conservative density sensitive 
limits throughout the observed region (SK89, logA^e.Sn ~ 2.45 
and logA^e.Oii ^ 2.4). Neglecting the lower limit the three in- 
nermost arcs, at 66"5 < r < 68"5, give an average value A^e.Sn = 
30 +4 cm"^, where the error is the standard deviation of the three 
values. Unlike for the other objects of this study we did not de- 
rive a density structure using the Abel transform. That procedure 
relies on a complete ionization in the integrated area; because of 
the low luminosity of the CS that does not appear to be granted. 

Sabbadin et al. (1985) deduce a mean density for the neb- 
ula of (A^e,Oii>^e,Sii) ~ 100 cm"-'. SK89 (using observational 
data from Osterbrock 1960 and Boeshaar 1974) likewise find 
A^e.on - 80 cmr^ and A^e.s,, = 170 cm^l Cuesta & Phillips (2000) 
present a density map using narrow -band imaging of [Sii]. In 
that map a distinct dichotomy is found in a relatively flat density 
gradient between the west and east sides. The average of their 
map, (A^e.Sii) = 590cm"-', appears unrealistic in view of our new 
- and previous - results. Kwitter & Henry (2001), finally, find 
A?e,s„ = 100+10cm-^ 



5.5. NGC 6826 - PN G083. 5+12.7 

While the central regions of NGC 6826 are often studied in the 
literature, this is not the case with the halo. Using the PPAK IFU 
we observed the full radial extent of the halo with two pointings 
('"halo - 65"). Compared with the object image in Fig. le (i.e., 
fig. 13 in CSSP03, also see GVM98 and Hajian et al. 2007) it 
is seen that the halo of this object deviates only slightly from a 
spherical shape, justifying our use of concentric arcs to bin spax- 
els; the two maps that we used are shown in Fig. 5. Note that the 
large diameter of the PPAK fibers prevents us from resolving the 
central region inside of the rim, which is why we use the radius 
given by CSSP03, - 6". The geometrical properties of the 
IFU also make it difficult to get perfectly overlapping arcs of 
spaxels in the two spatially offset tiles. Compare, for example, 
the spaxel organization of arc e in the two tiles. Intensities mea- 
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sured in the same arc, but in different tiles, could consequently, 
depending on the object properties, differ appreciably. Finally, 
because of the fiber separation and size, and since the obser- 
vations were carried out at a low airmass, we did not consider 
effects of DAR. 

To mention a few references, Kingsburgh & Armour (2001) 
measure an extinction coefficient of c = 0.25 4" NE of the CS. 
Barker (1988) finds c = 0.16 (3"N, 3"W), and KH98 measure 
c - 0.00-0.36 at three different central positions. In the central 
part of the nebula, we found a mean value of the extinction coef- 
ficient c = 0.24±0.11. Because of weak H-y emission it could not 
be well determined for the region outside the shell. Six distance 
estimates of NGC6826 are, moreover, reported by: MCW89 
(d = 2.3 kpc), Mendez et al. (1992, d = 1.9 kpc), CSSP03 
(d = 1.8 kpc), Pauldrach et al. (2004, d = 3.18 kpc), Phillips 
(2004, d=l.2 kpc), and Kudritzki et al. (2006, d = 2.6 kpc). We 
used J = 2.5 kpc (see Sect. 6.2). An effective temperature is given 
by, e.g., Kudritzki et al. (1997, T^g = 50 000 K) and Kudritzki 
etal. (2006, reff = 46000K). 

5.5.1. Surface brightnesses 

Radial surface brightness structures are presented in Fig. 19 
for seven emission lines. Note that the full extent of the limb- 
brightened halo is covered in H/?, Hy, and [Om] /14959. A com- 
parison of the intensities with the observed radial profiles of 
CSSP03 shows a good quantitative agreement (for HQr-i-[N n] and 
[O III] - diffusive light) - the amount of scattered Ught in the data 
is small. 

Outside the "halo boundary", for radii r>70", the three out- 
ermost measurements of Hfi are underestimated by about 10 per 
cent due to a varying continuum level around this einission line. 
Nevertheless, these data points are shown in order to demon- 
strate that emission is present. Finally, the helium and argon 
emission lines He 1/14471 and [Ariv] AA4111, 4740 are only de- 
tected in the inner part of the halo, out to a radius of about 

5.5.2. Electron temperature 

We measured an electron temperature out to a distance of r = 33" 
from the CS, using the radially binned intensities of [O m] /14363 
and [Oin] /t4959 - thus covering the inner halo. We present the 
derived radial temperature structure in Fig. 20. 

A steep temperature gradient is seen outwards of the shell, 
for r ^ 10", where the temperature increases from = 9200 K 
to 15 000 K within Ar ^ 12". The maximum is followed by 
a decreasing temperature in the three outermost arcs, where 
Tg > 12 000 K. Two measurements of the halo temperature 
are reported by Middlemass et al. (1989b, who measure - 
13 000:1:^000 Katr^ 46", hereafter MCW89) and Manchado & 
Pottasch (1989, who measure = 10600 ± 500K at r 36", 
hereafter MP89) - both values are indicated in Fig. 20. These 
values are measured farther out in the halo than what was pos- 
sible with the data presented here. The value of MP89 is about 
2000 K lower than we measured at r = 3 3 " . 

JQA87 present a temperature (and density) map of the cen- 
tral nebula. The agreement with our data is good. In addition to 
the halo values MCW89 and MP89 measure Te = 10400±750K 
and Te - 9 100 ± 100 K in the core, respectively. Both these pa- 
pers report on the apparent increasing gradient using the two 
respective values of the central parts and the halo. Additional 
values are reported by McKenna et al. (1996, who use data of 
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Fig. 20. Radial temperature structure of NGC6826, cf. Sect. 5.5.2. In 
addition the two halo measurements of MCW89 (A) and MP89 (❖) are 
indicated. For further details of figure properties see Fig. 9. 



Kaler 1986, Te = 11 llOK), KH98 (T^ = 8800-9000 ± 500 K), 
and Zhang et al. (2004, Te = 9350 ± 200 K). 

5.5.3. Electron density 

To mention a few references JQA87 provide a density map of 
the central nebula, using the sulfur line doublet, A^e.Sn- Using the 
same fines Barker (1988) and KH98 calculate a density at six and 
three different positions, respectively. MP89 determine average 
values from the Hy6 flux, A^e(core) - 2000 cm"^ and A^e(halo) = 
20cm-l Wang et al. (2004) calculate logA^e.On = 3.27loo^, 
logA?e,s„ = 3.28:0:06, logAre,cin, = 3.12:0-08, and logAre,Ar„ = 
3.33^0 2^ cm"-'. Zhang et al. (2004) find a value derived from 
spectra near the Balmerjump, logA'e,Bai = 3.4 + 0.1 cm"-'. We ob- 
served the two density sensitive argon lines, but in this case the 
heUum fine He 1/14713 is of similar strength as [Ariv]/I4711, 
which is why we could not deblend the two fines accurately. 
Note that the expected electron density also is below the con- 
servative sensitivity limit of the argon fine doublet (SK89, 

l0gAfe,Ariv^3.3). 

Following the approach of PS90 for NGC 6826 we also used 
an Abel transform on the radial structure of H/3 in order to calcu- 
late a density structure. PS90 fix the density to A'e = 1000 cm"^ at 
the radius r = 12". We derived a similar structure under the addi- 
tional assumptions of a distance of J = 2.5 kpc and a temperature 
of Te = 12 000 K. With these values we also get A'e('- = 12") = 
1000 cm"^, which at this radius is well within the azimuthal den- 
sity range reported by JQA87 (800 < A^e.Sn ^ 2000 cm"^). The 
result is shown in Fig. 21, also see Sect. 6.2 for a discussion on 
the derived mass loss rate. 



6. Discussion 

In this section we first discuss the appearance of hot halos in 
Sect. 6.1. Thereafter we derive mass loss rates and core and halo 
masses for three of the objects of this study in Sect. 6.2. 
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Fig. 19. Radial surface brigtitness structures of NGC6826, cf. Sect. 5.5.1. From the top the panels show seven emission line intensities of four 
elements: a) hydrogen, b) helium and argon, and c) oxygen. For further details of figure properties see Fig. 8. Data of lines drawn in gray were 
measured in tile 1, and those in black in tile 2 (see Table 1). The vertical gray bars indicate the location of the rim at = 6", and the shell-halo 
transition at rjheu = 12". Note the vertical offset between black and gray lines in the inner halo, for 15" <r< 25". This offset is due to deviations 
from synmietry. 



6. 1. Occurrence and understanding of hot halos 

A knowledge of the electron temperature structure is important 
when determining the ionization structure and abundances of el- 
ements in PNe, including the halo. In our halo observation data 
we detected emission lines of hydrogen, oxygen, neon, and oc- 
casionally also helium and argon in the inner halo. Halo tem- 
peratures can unlikely be determined by any other method than 
through the oxygen [O m] emission lines. 

MCW89 (NGC6826), MP89 (NGC6826 and NGC6543, 
also see Martin-Gordon et al. 2007, for the latter), MCWH91 
(NGC 7662), and Monreal-Ibero et al. (2005, NGC 3242) present 
temperatures calculated at one position in the core and one in 
the halo for the respective object. Radial temperature structures 
which seem to indicate high halo temperatures for NGC 1535 
and NGC 3242 are presented by Krabbe & Copetti (2005, see 



figs. 5 and 13 therein). Compared to these studies we derived 
densely sampled temperature structures for four out of five neb- 
ulae, reaching from the FN core and out into the halo. In three of 
the four temperature structures we found a hot halo: NGC 7662 
(Fig. 9), NGC 6826 (Fig. 20), and IC 3568 (Fig. 14). The figures 
show that the region is narrow, where the temperature increases 
from the cooler PN core to the hotter halo. Note that de-projected 
(real) temperature structures are expected to be different, and 
temperatures can be even higher. Very high temperatures could 
be important, as Ha becomes collisionally excited when temper- 
atures are higher than 20 000 K (Stasinska 2002). 

Early studies, which are based on hydrostatic model formu- 
lations, explain hot halos as a consequence of radiation hard- 
ening or shock heating mechanisms (see, for example, MP89; 
MCW89; MCWH91). Due to high densities a PN core is in gen- 
eral in thermal equilibrium, but the tenuous halo is not. Non- 
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Fig. 21. Radial density structure of NGC 6826, cf. Sect. 5.5.3. The thick 
gray line shows a fit of the density to the power law A^e " for the 
region 15" <r<33". For further details of figure properties see Fig. 11. 



equilibrium models, which account for hydrodynamic effects, 
are first suggested to explain high temperatures in the halo by 
Marten (1993, also see Marten 1995 and Mellema & Frank 1995; 
the models are based on the formalism introduced by Marten 
1994 and Marten & Szczerba 1997) and Tylenda (2003). With 
such models the ionization front, which is responsible for the 
high temperature in the halo, moves through the envelope in 
a short time period (about lO^-lO"* years long). These models 
suggest that high temperatures in the halo are a fairly common, 
although transient, phenomenon. 

In Figs. 22a-c we show a sequence of radial PN structures 
that were calculated using a time-dependent radiation hydrody- 
namic model, where T^ff ^ 80 000 K and L ^ 6000 Lq. The start- 
ing model, for the calculations of the hydrodynamic evolution 
through the PN stage, is taken from Steffen et al. (1998, see 
fig. 21 therein). We selected this model since it shows how the 
R-type ionization front moves through the former AGB wind, 
thereby giving rise to an observable halo. We also wanted to 
show qualitatively how the high temperature occurs during the 
ionization of the halo - the model is not a fit to the observa- 
tion data. An R-type ionization front is moving outwards and 
ionizes the halo, see the blue solid line that shows the electron 
density. In only 125 years the spectroscopic"^ electron tempera- 
ture in the central halo (where 7 S r 5 12 x lO'^cm) changes 
from Te ^ 26 000 K, in Fig. 22a, to ^ 19 000 K, in Fig. 22c, 
see the red line marked with asterisks. Note that in compari- 
son to the spectroscopic temperature, which is projected onto 
the line of sight, the real temperature (dashed line) can be about 
lO'^K larger in the same region. The spectroscopic temperature 
in the halo is strongly affected by the cool and dense wind in- 




* The spectroscopic temperature is derived from the computed line 
profiles of [O m] as a function of the distance from the CS. 
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Fig. 22. A radiation hydrodynamic model sequence showing the evo- 
lution of a radial PN structure as the halo is being ionized. Between 
panels a and c the model ages by 125 years, and the effective temper- 
ature and the luminosity of the CS change from T^g.^ = 76 800 K and 
La = 5850 Lq to Teff.^ = 80 200 K and L, = 5790 Lq. The four lines 
show: the total number density of atomic nuclei (black solid line), the 
electron density (blue solid line), the real electron temperature (dashed 
line), and the spectroscopic temperature (red line connected with aster- 
isks). The density ordinates are logarithmic (left hand side axes). The 
halo is the region between the shell-halo transition, at rsden 3 x 10' ' cm 
andrhi,io-12x 10" cm. The radial domain 12 x 10"<r<18x 10" cm 
is the wind interaction region of the last thermal pulse. Mass loss occur- 
ring before the last thermal pulse is found where r> 18 x 10" cm. 



teraction region at 12 ^ r ^ 18 x lO'^ cm, causing the rapid de- 
crease. A use of the spectroscopic temperature to calculate abun- 
dances would therefore give incoiTect values. Figure 23 shows 
the surface brightnesses of Hfi and [O iii] /15007. Note the inten- 
sity maximum in the halo for 13Sr$17x 10'^ cm, that is also 
visible in Fig. 19d for NGC 6826. 

6.2. Deriving mass loss rates for winds at fhe tip of tlie AGB 

Starting with an electron density (Ng), that we derived using an 
inverse Abel transform, we first calculate a mass density (p). 
With a mass density structure and a velocity structure it is then 
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Fig. 23. The surface brightness of HyS (red dashed line) and [O m] /15007 
(solid line) calculated using the physical structure shown in Fig. 22c. 
For further details see Sect. 6.1. 

straightforward to calculate a radially dependent mass loss rate 
(M). Under the three assumptions of: full ionization, a simple 
chemistry of 90% H and 10% (singly ionized) He, and a velocity 
structure u(r) we get, 

p (r) = [0.9 X 1 + 0. 1 x4]mMr) and M (r) = Anrp (r) u (r) , 

where mu is the atomic mass constant. The velocity structure 
is in general not known, which is why we instead provide the 
following property, 

M (r) /« (r) = Anp-p (r) , 

that can be multiplied by any radially dependent velocity to get 
the mass loss rate. In this study we assumed in every case a con- 
stant velocity, u (r) = 10 km s"' . We present radial mass loss rate 
structures - that in every figure are indicated by the line that is 
connected by square symbols - for NGC 7662 (Fig. 1 1), IC 3568 
(Fig. 15), M2-2 (Fig. 18), and NGC 6826 (Fig. 21). We did not 
consider NGC 3587 since it is an old object that is not expected 
to contain any information of the AGB wind anymore. Note that 
densities and masses depend on the assumed distance (d) to the 
object: mass loss rates scale with the distance as M oc d^^^ and 
integrated masses as Mocd^^^. 

Only parts of a halo contain information of the mass loss 
rate on the AGB. In the region that is close to the shell pho- 
toionization processes have affected the wind structure, and in 
the outer parts of the halo the wind structure has been affected by 
interactions between the wind and previous mass loss episodes. 
For three out of the four objects, where we measured density 
structures, we observed only the inner halo. The resulting "trun- 
cated" intensity structure of these objects introduces an error 
into the derived density, that is above a few per cent only in 
the outermost parts. In the central nebula (i.e., the PN core) 
an asymmetric structure (NGC 7662) or a low spatial resolution 
(NGC 6826) makes the density inaccurate there, which is why 
our core masses likely are uncertain for these objects. 

We present the results of our study in Table 2. In every case, 
but M2-2 (where there is no comparison), we based the dis- 
tances on the Shklovsky distances (and luminosities) given by 
Gorny (priv. comm., see CSSP03). The four PNe are in a simi- 
lar evolutionary state and should have similar luminosities, and 
we therefore scaled the distances to get the same luminosity, 
L ^ 6000 L©, that corresponds to a CSPN mass of ^ 0.6 Mq. 
When we calculated halo masses we used the entire radial extent 



of the density structure we had outside the shell. We only ob- 
served the full radial extent of the halo for NGC 6826, the halo 
masses of the remaining objects are lower limits only. When we 
calculated core masses we used the region between the rim and 
the shell. Our distances are larger than those given in the liter- 
ature, and our masses are consequently also larger. In Table 2 
we also compare our core and halo masses with literature values 
(Mcore.ref and Mhaio.ref)- Sincc adopted distances differ between 
references we scaled the respective value to use our distances. 

Our masses should be in best agreement with PS90, since we 
used the same method as they do. In comparison Hippelein et al. 
(1985) assume the density dependence n oc r"^, and the other 
three papers use a constant density in the halo. If we truncate 
the density structure at r = 70", that is the outer boundary used 
in PS90, we get the halo mass Mhaio - 1.5 M©, a value that is 
in good agreement with Mhaio,ps90 = 1.3 Mq. Compared to PS90 
we, furthermore, did not use the entire halo surface to determine 
the density structure, and the spatial resolution of our data is 
lower, in particular in the PN core. We attribute the remaining 
deviation in the outcome to these differences. 

A power law fit of the density with the radius (p oc r^") 
should, if the mass loss rate and outflow velocity are constant, 
give the exponent a = 2. For each object we measured several 
values of a in the same region where we calculated a mass loss 
rate. The upper and lower limits of these values are presented 
in Col. 5 of Table 2. In every case a is greater than 2, and in- 
creases to Q- ^ 5-8 in the inner part of the halo - indicating a 
marked increase of mass loss during the final outflow stage on 
the AGB. Fits of a for separate regions are shown in the density 
plots of the respective object; the value of NGC 6826, a - 4.2, 
agrees well with that derived by PS90, a = 4. Note, however, that 
the ionized AGB wind is accelerated by its own pressure gradi- 
ent and is no longer exactly described by the original power law 
index a' of the neutral wind (see Schonberner et al. 2005). 

The property ^, that is the ratio between the minimum and 
maximum measured mass loss rates, varies between 3.6-6.7. 
Note that both a and ^ are independent of the distance. Despite 
the relatively low luminosity (of only ~ 6000 Lq) and increas- 
ing effective temperature of the PPNe moving off the tip of the 
AGB, the maximum mass loss rates average at M^ax - lO"'* Mq. 
Because of the presence of scattered light, in that part of the halo 
which lies the closest to the shell, mass loss rates are there over- 
estimated by up to about 30 per cent (Sect. 4. 1); hence real slopes 
are less steep. Moreover, a use of shorter distances, which are 
smaller by a factor of 2, result in a factor 3 lower average mass 
loss rate. The four mass loss rates are also similar, within a factor 
of 2. These mass loss rates are important constraints on the the- 
ory of wind formation, which so far does not treat winds on the 
tip of the AGB satisfactorily. Results of earlier studies which in- 
dicate an increasing mass loss rate are reported by MP89, PS90, 
and Schonberner et al. (2005), also see Kwok et al. (2002). 



7. Conclusions 

Surface brightnesses of emission lines in the halo of a PN can 
be several thousand times weaker than in the PN core. Therefore 
it is a demanding task to derive physical properties of the halo. 
Previous studies of PNe, which focus on the halo, rarely derive 
physical properties at more than one position. For this paper we 
wanted to make an observational study of PNe in the galactic 
disk, where we put the emphasis on the halo and the shell-halo 
transition regions. Our goal was to derive densely sampled tem- 
perature structures extending from the PN core and as far out in 
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Table 2. Calculated mass loss rates and core and halo masses, cf. Sect. 6.2. Object names and distances (d) are given in the first two columns. 
Columns 3-8 specify six properties: the effective temperature (T^g), the radial range Ar = rmm-''max where we calculated a mass loss rate, the range 
of exponents a (of fits to oc r'") in the radial range Ar where OTmin was calculated just inside r^ax and amax just outside r^iin, the maximum mass 
loss rate in the same range (Mmax), the maximum-to-mirrimum mass loss rate ratio in the same range (^), and the resulting mass of the core and 
the halo in the region where we calculated a density (Mcore, Miaio)- The halo mass is printed in italics where the full radial extent of the halo was 
not covered. The subscript index denotes the age of the radial extent where we integrated a halo mass (specified in thousands of years). A constant 
outflow velocity of «= lOkms"' was assumed throughout all objects. The comma separated core and halo masses in Col. 9 (Mcore_ref> A^haioref) are 
based on the values given in the respective reference, that are specified in Col. 11. The original masses, that used the distances specified in Col. 10 
(rf,rf), were scaled to use our distances (d, Col. 2). References: (1) Hippclcin ct al. (1985). (2) MCW89, (3) MP89, (4) PS90, and (5) MCWH91. 



Object 


d 




mill ^ max 




10^M„„ 




Afcore • -^halo 


-^corejef • -^halo,ref 


disf 


Ref. 




[kpc] 


[io^k] 


[arcsec] 




[Moyr-'] 




[Mq] 


[Mq] 


[kpc] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC7662 


1.7 


100 


15-36 


6.6-2.4 


10 


6.7 


<0.64 : 0.55,7 


0.053:0.32, 0.19:0.26 


1.15, 1.5 


(1), (5) 


IC3568 


4.3 


55 


16-33 


7.9-2.4 


7.4 


3.6 


0.77 : 0.9h^ 








M2-2 


4.0 


80 


10-14 


7.3-2.6 


12 


4.7 


0.50 : 0.8825 








NGC6826 


2.5 


46 


6.8-13 


5.5-3.4 


7.1 


5.5 


0.50 : 1.976 


0.56:1.8, 0.52:1.2, 


1.02, 2.3, 


(1), (2), 


















0.59:1.1,0.50:1.3 


1.0, 1.54 


(3), (4) 



the halo as possible. And also to study the final mass loss history 
on the AGB. 

We carried out the observations using the IFS instrument 
PMAS. Through this approach we could account for several is- 
sues - in our data analysis - which treatments are critical to 
the outcome when evaluating steep intensity gradients, such as 
found in PNe. Next we summarize our results of four such is- 
sues. First, in every case, but one (where it was not necessary), 
we corrected individual spectra in the reduced data for DAR. 
Depending on observing conditions we thereby avoided errors in 
derived temperatures of about a thousand Kelvin. Second, in or- 
der to save valuable observing time we developed a new method 
to subtract telluric emission lines from the spectra, instead of 
using separate sky exposures. A telluric emission line of key im- 
portance to our studies is the mercury line Hg/14358 that often 
is strong and lies adjacent to the temperature sensitive, and often 
very weak, oxygen line [O iii] /14363. At the rather low spectral 
resolution of our observations a careless subtraction of Hg /14358 
may lead to temperatures which are several thousand Kelvin too 
large. Our novel method, where telluric lines and line wings, are 
fitted with a set of fines of fixed proportions - together with the 
object lines - is useful and general for as long as the telluric and 
the object fines are resolved. 

Third, in order to increase the signal-to-noise of individual 
emission lines we binned spectra of individual spatial elements 
on the IFUs along concentric arcs - under the assumption that 
the halo structure is spherically synametric. The arcs were cen- 
tered on the CS. With this method we extracted intensities at in- 
creasing radii - at a spatial separation and width of Ar > 0'.'5-3'.'6, 
depending on the choice of the IFU and line strengths. Fourth, 
and finally, we made a study of the contribution of scattered light 
to the data when using PMAS and the 3.5m telescope at Calar 
Alto. The results showed that some scattered light is present 
in the inner halo. Its contribution at any radial point is moder- 
ate, ^ 30%, and mainly comes from bright nearby regions; from 
within a distance of about 4" when using the I'.'O IFU. Our ob- 
servations of faint halos consequently show real structures. 

For the extended and old nebula NGC 3587 we could deter- 
mine one value of the electron temperature in the outer parts of 
the nebula that agrees with earlier values in the literature. Of the 
remaining four objects in our study, we found that M 2-2 has a 
halo radius that is more than twice as large as reported earlier 
(''halo - 20"). For IC 3568 we found a halo (rhaio ^ 25") where 
this had not been detected earfier. We, moreover, calculated ra- 



dial temperature structures covering the PN core and parts of 
the halo for NGC 7662, IC3568, M2-2, and NGC 6826. While 
the structure of M 2-2 showed a large temperature in the cen- 
ter, and an otherwise constant value, the remaining three struc- 
tures showed a steep radial temperature gradient with sigiufi- 
cantly higher values in the halo. Hot halos are described with 
the help of hydrodynamic models as transients, which in a short 
time move through the halo. Since matter in halos is rarefied it 
takes a long time before a equiUbrium conditions are established. 

We derived electron densities in the shell for NGC 7662 
(A'e.Ariv) and M 2-2 (A'e.s n)- Using radial surface brightness struc- 
tures, and assuming a constant radial velocity, we additionally 
derived density and mass loss rate structures, and masses, for 
NGC 6826, NGC 7662, IC3568, and M 2-2. AU mass loss rate 
structures show an increasing outflow in the final stage of mass 
loss on the AGB. Comparing the four ratios of maximum to min- 
imum mass loss rates these are about 4-7, and within a factor of 
two the maximum mass loss rate is Mmax - 10 Moyr"' . A re- 
markable value in view of the relatively low luminosity of these 
objects that is close to 6000 L0. It would be of great interest to 
see what the outcome would be of theoretical models of mass 
loss, which assume the physical properties at the tip of the AGB. 

FinaUy, a logical foUow-up of this study is to calculate abun- 
dances of the AGB wind in the inner halo. This is, however, only 
possible for a few elements - such as helium, oxygen, neon, and 
maybe also argon - since intensities are so weak. The results 
of our study should, moreover, be valuable to comparisons with 
objects in the halo and the Magellanic Clouds. In a new study 
we are working on the detection and measurement of halos in a 
sample of PNe in the Small and Large MCs, in order to study 
closer the influence of the metaUcity on the physical structure. 
The spatial resolution, which can be achieved with these objects 
is naturally lower than with galactic disk objects, since the two 
MCs are more distant. 

Acknowledgements. CS. is supported by DFG under grant number SCHO 
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Fig. 24. Signal-to-noise maps of a mosaic of 
tlie three observed tiles of NGC 7662, for two 
emission lines. In comparison to Fig. 7 this fig- 
ure shows the inverse relative surface bright- 
ness error for the same tiles. Values are log- 
arithmic. From the top the two panels show: 
a) [O m] ^5007 and b) [O in] /14363. For further 
details see Sect. 5.1.1. 
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Fig. 25. Binning maps used in the analysis of the observations of IC 3568, cf. Sect. 5.2. The map on the left hand side was used with emission lines 
in the blue part of the spectrum, for wavelengths shorter than /l^f = 5 100 A. The map on the right hand side was similarly used with emission lines 
of wavelengths longer than /l,ef. Each spaxel (square) represents an area of I'.'O x I'.'O on the sky. For further details of figure properties see Fig. 4 
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Fig. 26. Binning maps used in the analysis of the observations of M2-2, cf. Sect. 5.3. Each spaxel (square) represents an area of ff.'S x O.'S on 
the sky, and data is here present out to a distance of 20" west of the CS (•). The map in the top panel is used with all emission lines bluewards 
of /Iref = 5050 A, and the map in the bottom panel with all emission lines redwards of this wavelength. The data in the red wavelength range was 
observed in a subsequent night, at a slightly different angle, which is why it requires a different pattern mask with the DAR correction. For further 
details of figure properties see Fig. 4 
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Fig. 27. Surface brightness maps (left) with corresponding signal-to-noise maps (right) of tiles 1 & 2 of M2-2, cf. Sect. 5.3.1. From the top the 
panels show four emission lines of decreasing wavelength: a) [Oiii] /15007, b) H/?, c) [Oiii] /i4363, and d) [Neiii] /13869. For further details of 
figure properties see Fig. 7. Ansae are partly visible at (RA,Dec)-offsets (-3'.' 5,-1'.' 5) and (-5,-r.'5). Note the low signal-to-noise of the temperature 
sensitive oxygen line [O iii] /i4363 in panel c, already at a distance of r ^ 4" from the CS (■A'). 
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Fig. 28. Binning map used in the analysis of the observations of NGC3587 for the red wavelength range, cf. Sect. 5.4. Each spaxel (square) 
represents an area of I'.'Ox I'.'O on the sky. The CS is located 66" to the west (right) of tile 1. For further details of figure properties see Fig. 4 
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Fig. 29. Radial surface brightness structures of NGC3587, cf. Sect. 5.4.1. From the top the panels show intensities of 16 emission Unes of six 
elements: a) hydrogen, b) helium and neon, c) oxygen, d) sulfur and oxygen, and e) nitrogen. The locations of the rim at = 85", and the 
shell-halo transition at rsheii = 105", are indicated with gray vertical lines. Data of lines drawn in black were observed in April 2006, data of Unes 
drawn in gray in February 2004. For further details of figure properties see Fig. 8. 



